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PREFACE 


The  Defense  Nuclear  Agency  has  collected  and  printed  the  attached  papers 
from  the  February  25-27  1986  Global  Effects  review  as  a  service  to  the 
community.  The  Defense  Nuclear  Agency  takes  this  opportunity  to  express  Its 
gratitude  to  the  numerous  participants  In  the  Global  Effects  review. 

The  technical  papers  enclosed  Include  all  those  which  were  received  by  DNA 
prior  to  the  closing  date  of  28  April  1986.  Where  papers  are  missing  their 
place  Is  occupied  by  the  abstract  received  prior  to  the  meeting. 

The  Inclusion  of  a  paper  In  this  proceeding  does  not  necessarily  Imply 
endorsement  of  the  results  of  the  research  reported  or  conclusions  which  might 
be  drawn  from  that  research.  It  Is  the  opinion  of  the  Defense  Nuclear  Agency 
that,  while  good  progress  Is  being  made  In  Improving  our  understanding  of 
Global  Effects,  the  results  to  date  are  tentative  and  preliminary  and  should 
not  be  used  for  planning  beyond  the  planning  of  future  research. 
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GLOBAL  MODELING  &  APPLICATIONS 
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Earth  and  Space  Sciences  Division 
Los  Alamos  National  Laboratory 


PREVIOUS  VALIDATION  STUDIES 

AEROSOL  RESIDENCE  TIMES 

-  lower  and  upper  troposphere  (quantitative) 
[where  smoke  would  be  injected] 

-  lower  stratosphere  (qualitative) 

VERTICAL  TRANSPORT  OF  AEROSOLS 
(qualitative) 

-  rapid  vertical  mixing  within  troposphere 

-  slow  penetration  from  troposphere  into 
stratosphere 

-  slow  vertical  mixing  within  stratosphere 


QUESTIONS  ARISING  FROM 
NUCLEAR  WINTER  SIMULATIONS 


Vertical  transport  of  aerosols  due  to  atmospheric 
heating  by  aerosol  absorption  of  solar  radiation 
no  suitable  analogue 

Poleward  migration  of  aerosols  injected  in 
midlatitudes 

Arctic  haze ,  heavy  methane  tracer ,  volcanic  and 
atmospheric  weapons  test  debris 

Residence  time  of  aerosols  in  stably  stratified  layers 
volcanic  and  weapons  test  debris  (stratosphere) 
Arctic  haze  (polar  inversion  layer) 

Rate  of  dispersion;  inhomogeneity;  "patchiness" 
heavy  methane  tracer  (troposphere) 
volcanic  and  weapons  test  debris  (stratosphere) 


...  .V 


O 

z  co 
<  c c 

Offi 

HI 

PQ 


CO 
«  LU 

OQ 


Q  -J 

DC  3 


>z 

>< 

S5f 


O) 

> 

0 

A 

0 

■  Hi 

i_ 

B^B 

C 

0 

n 

X 

CO 

CO 

o 

2 

LL 

2 

CO 

> 

5 

CO 

0 

CO  5 
o 


Ui  _ 
c  TO 

v  ° 
.2  O 

I  o 


a 

co 


3: 

£  ° 
0 


a 

co 

o 

a 

o  *- 

h| 

o 


0 

0 

E 

o 

o 

0 

c 

> 

0 

o 

0 

o 

o 

'5 

c 

0 

E 

3 

0 

0 

O 

DC 

1- 

DC 

CO 

c 

a> 

i  w 

o  __ 


Ol  .2 


a  co 

co  a 
co  co 


servability  Visual  +  Sampling  Sampling  Remote  + 

Sampling  only  only  Sampling 


ANALYSIS  OF  OBSERVATIONAL  DATA 
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EULER1AN  TRANSPORT  MODEL 


Calculates  tracer  concentration  at  discrete  grid 
points 

Has  same  resolution  as  underlying  GCM 
Limitation:  number  of  grid  points  affordable  in  GCM 


LAGRANGIAN  TRACER  MODEL 


Tracer  is  represented  by  a  large  number  of  discrete 
particles 

Particle  coordinates  not  restricted  to  GCM  grid 
points 

Limitation:  number  of  particles  needed  to 
adequately  represent  local  concentrations  on  grid 


ARCTIC  HAZE 

Composed  of  industrial  pollutants  (heavy  metals, 
sulfates,  organics,  gases) 

Major  source  region:  Eurasia,  especially  USSR 

Pollutants  from  eastern  U.S.  are  washed  out  in 
storm  track  over  North  Atlantic 

Pollutants  transported  directly  to  Arctic  via  lower 
troposphere 

Haze  found  from  surface  up  to  at  least  6  km 
Maximum  concentration  in  winter  to  early  spring 

Reasons  for  winter  maximum: 

-  pollutants  accumulate  in  stable  Eurasian  air 

-  little  precipitation  in  Arctic  in  winter 

-  vigorous  circulation  enhances  transport 


Model  correctly  portrays  (qualitatively) 

•  accumulation  of  pollutants  over  Eurasia 

•  poleward  migration 

•  scavenging  of  pollutants  from  North  America 

•  seasonal  contrast  in  haze  concentration 


AEROSOL  RESIDENCE  TIMES 


TABLE  8-5  f  \\ n ^ \ 

Residence  lime  of  atmospheric  aerosol  particles  at  various  levels  in  the  at 

mosphere. 


LOWER  TROPOSPHERE  <1  WEEK 
UPPER  TROPOSPHERE  -  1 .5  WEEKS 
LOWER  STRATOSPHERE  -  3  MONTHS 


ilometers 


%of  total  activity 
in  5  kft  increments 

6/27/73  _  8 

Test  Assumed 

-  I  -  2  - 


June  27,1973 


Assumed 


Estimated  tropopause  height 
of  the  6/27/73  test 


PERCENT  OF  TOTAL  ACTIVITY 

Fig.  8.  Comparison  of  the  assumed  initial  vertical 
activity  distribution  for  a  3  MT  total  yield  nuclear 
test  in  temperate  latitudes  (based  on  the  4  earlier 
high  yield  Chinese  tests)  with  the  estimated  initial 
vertical  activity  distribution  of  the  June  27,  1973 
Chinese  test  of  2  to  3  MT  total  yield. 
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Sedlacek  et  al.:  Stratospheric  SO«“  and  Volcanic  Activity 


INTERCEPT  =  440  PPBM 


K 


1/e  FOLDING  TIME  -  11.2  MONTHS  FOR  VOLCAN  FUEGO 

OH - -I - 

0  05 

YEARS  AFTER  ERUPTION 

Fig.  35.  Average  ‘stratospheric’  sulfate  concentration  minus  preceding  ‘ambient’  concentration  versus  time  from  erup¬ 
tion  to  middate  of  sampling  period.  The  least  squares  fitted  e-fold  removal  rate  for  Fuego  1974  is  1 1.2  months.  Only  the  first 
five  data  points,  indicated  by  circles,  were  used  to  determine  the  least  squares  fit  due  to  the  low  volcanic  sulfate  to  ‘ambient' 
sulfate  ratio  for  the  last  two  points.  For  Mount  St.  Helens  the  contribution  of  Sierra  Negra  preceding  and  subsequent  Mount 
St.  Helens,  Gareloi,  Alaid,  and  Pagan  eruptions  make  an  e-fold  removal  rate  calculation  meaningless. 
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STRATOSPHERIC  RESIDENCE  TIME 


Some  disagreement  among  observations,  but 
recent  data  indicates  1 /e-times  of  6-12  months  in 
lower  stratosphere 

Observations  show  comparable  residence  times  in 
stratosphere  for 

-  gases  (HTO) 

-  radionuclides 

-  volcanic  sulfates 

Radionuclide  data  suggests  no  removal  for  a  long 
period  (3-6  months)  following  lower  stratospheric 
injections 

Eulerian  transport  model  gives  1/e-times  in  lower 
stratosphere  of  2-3  months 


Three  releases  of  1  kg  CD4  in  midtroposphere  near 
55S,  165E  (January,  June,  October  1984) 

Precisely  characterized  source:  time,  place, 
altitude,  amount.  Inert  tracer. 

Samples  collected 

-  on  airplane  flights  between  Antarctic  stations, 
and  to/from  New  Zealand  (1  hour  averages) 

-  at  surface  stations  (3  day  averages) 

-  for  60  days  after  release 

Detection  limit  -  10“17  g  CD4  /  g  air 

Samples  above  limit  found  out  to  20  days 
(longest  data  series  for  a  tropospheric  tracer?) 

Principal  shortcoming:  limited  sampling  precludes 
complete  picture  of  large-scale  distribution. 
Modeling  comparisons  may  help. 


MODEL  VS  OBSERVATIONS 


Eulerian  transport  model 

-  gives  reasonable  qualitative  tracer  behavior 

-  but  concentration  is  diluted  too  rapidly 


•  Inherent  inability  of  Eulerian  model  to  simulate 
"point  source" 


•  Dispersed  source  at  t=0  may  make  simulated 
tracer  prematurely  susceptible  to  dispersal  by 
large-scale  eddies 

•  Excessive  lateral  diffusion  of  Eulerian  tracer 
distribution  stretched  by  "geostrophic  turbulence 


GEOSTROPH1C  TURBULENCE 

Nonlinear  large-scale  eddies  in  rotating  stratified 
fluid 

Essentially  two-dimensional  (D  «  L) 

Energy  transfer  to  larger  scales  causes  stretching 
and  deformation  of  tracer  distribution  without 
mixing  (unlike  small-scale  3-D  turbulence) 


time  to  time  t, 

Fig.  4.  -  Mixing  in  a  two-dimensional  phase  space. 
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GEOSTROPHIC  TURBULENCE 


Growth  to  a  rather  large  scale  is  required  before 
deformation  becomes  rapid  and  pronounced. 

However,  if  start  with  large  distributed  source, 
deformation  can  begin  immediately. 

For  point  source,  time  required  for  distribution  to 
reach  the  critical  scale  depends  on  small-scale 
processes  that  may  not  be  well  represented  by 
coarse  resolution  of  GCM. 


LAGRANGIAN  TRACER  MODEL 


Use  to  study 

•  Dispersal  from  "point  source"  release 

•  Relative  importance  for  large-scale  mixing  of 

-  particle  distribution  deformation  (stretching) 
by  geostrophic  turbulence 

-  small-scale  motions  (unresolved  by  GCM) 
that  cause  spreading  relative  to  "stretched" 
distribution 

-  shearing  of  distribution  (in  the  vertical)  by 
vertical  gradients  of  the  horizontal  winds 

-  initial  vertical  and  horizontal  extent  of  particle 
distribution  ("exposure"  to  gradients  in  wind) 

•  Long  stratospheric  residence  times 


CONCLUSIONS 


Eulerian  transport  model 

-  is  suitable  for  distributed  sources 

-  captures  basic  features  of  Arctic  haze 

-  simulated  poleward  migration  from 
midlatitudes  agrees  with  observations 

-  displays  deformation  of  tracer  distribution 
characteristic  of  geostrophic  turbulence 

-  but  coarse  resolution  causes  too  rapid 
dispersal  compared  to  point-source  data 

-  underestimates  aerosol  residence  time  in 
lower  stratosphere  (recent  measurements 
give  6-12  months  or  longer) 

•  Most  promising  analogues  for  quantitative 
transport  comparisons: 

-  heavy  methane  (troposphere) 

-  radionuclides  (stratosphere) 


The  Hydrological  Response  to  Summertime  Injections 
of  Smoke  into  the  Atmosphere 

Steven  J.  Ghan,  Michael  C.  MacCracken,  John  J.  Walton 
Lawrence  Livermore  National  Laboratory 
January  1986 

Studies  of  the  summertime  temperature  response  to  varying  amounts  of  smoke  indicate 
that  sharp  cooling  does  not  occur  until  sufficient  smoke  has  been  injected  to  convectively 
decouple  the  surface  and  troposphere.  For  smaller  amounts  of  smoke,  relatively  small 
temperature  changes  may  result,  and  even  modest  warming  may  occur  as  a  result  of  a 
reduced  planetary  albedo.  Although  convective  processes  may  still  be  tieing  the  surface 
and  troposphere  together  for  these  smaller  smoke  amounts,  the  rate  of  convection  drops  in 
direct  proportion  to  solar  absorption  by  the  smoke  particles.  Thus,  convection  is  suppressed 
for  small  smoke  loadings  as  well  as  large  ones.  Because  much  of  the  normal  summertime 
precipitation  is  convective  in  nature,  we  conclude  that  significant  drought-like  reductions 
in  land  precipitation  are  possible,  even  for  summertime  tropospheric  smoke  injections  small 
enough  to  cause  only  small  temperature  changes,  or  even  a  surface  warming.  Results  from 
experiments  with  a  tropospheric  general  circulation  model  under  a  variety  of  summertime 
smoke  injections  are  presented  in  support  of  this  conclusion. 


The  Fate  of  Smoke:  Sensitivity  Studies  with 
a  Global-Scale  Atmospheric  Model 

Steven  J.  Ghan,  Michael  C.  MacCracken,  John  J.  Walton 
Lawrence  Livermore  National  Laboratory 
January  1986 

The  relative  importance  of  several  smoke  removal  processes  has  been  considered  using 
a  tropospheric  general  circulation  model.  Smoke  particles  have  been  divided  into  two  cat¬ 
egories:  those  smaller  and  those  larger  than  one  micron  in  diameter.  Removal  processes 
for  each  aerosol  include  precipitation  scavenging  and  dry  deposition.  Small  particles  can 
coagulate  to  become  larger  particles.  Scavenging  is  found  to  be  the  most  efficient  removal 
mechanism  for  the  large  particles.  For  sub-micron  particles,  the  scavenging  rate  is  based  on 
removal  rates  inferred  for  Arctic  soot,  primarily  by  ice  processes.  For  these  small  particles, 
coagulation  is  most  important  at  high  concentrations,  while  scavenging  and  dry  deposition 
are  comparably  important  at  lower  concentrations,  depending  on  the  specified  parameter 
values.  The  lifetime  of  the  smaller  particles  is  significantly  longer  than  for  the  large  parti¬ 
cles,  both  because  of  less  efficient  scavenging  and  because  of  a  lower  dry  deposition  velocity. 
The  importance  of  accounting  for  sub-grid  scale  patchiness  of  precipitation  in  reducing  ef¬ 
fective  scavenging,  and  the  smoke-induced  stabilization  of  the  planetary  boundary  layer  in 
reducing  dry  deposition  are  assessed. 

This  work  was  performed  under  the  auspices  of  the  U.S.  Department  of  Energy  by  the 

Lawrence  Livermore  National  Laboratory  under  Contract  W-7405-Eng-48. 
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Figure  22.  Lifetimes  of  large  smoke  particles  foi 
deposition,  precipitation  scavenging,  and  coag 
particles  the  coagulation  lifetime  is  not  defined 
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MERIDIONAL  CIRCULATION 

APRIL  CONTROL  (Days  16-20) 


Smoke  + 


Physical  Processes  to  Incorporate  or  Improve 
In  Three-Dimensional  Models 
For  Nuclear  Winter  Research 


Radiative  Transfer 

■  Infrared  effects  of  smoke  and  dust 

■  Scattering  of  solar  radiation  by  smoke  and  dust 

■  Refined  optical  properties  of  smoke  and  dust 

(size  distribution,  refractive  indices,  non-sphericity) 

Surface  Temperature  Calculation 

■  Diurnal  cycle  of  heating 

•  Explicit  parameterization  of  soil  heat  storage 

•  Vegetation  effects 

Planetary  Boundary  Layer  Representation 

•  Stability  dependent  vertical  mixing 

■  Inclusion  of  low  level  clouds  and  fog 


Physical  Processes  to  Incorporate  or  Improve 
In  Three-Dimensional  Models 
For  Nuclear  Winter  Research 


Sub-grid  Scale  Vertical  Transport 

■  Internally  consistent  convective  transport  methods 
for  smoke  and  dust,  as  well  as  heat  and  moisture 

•  Stability  dependent  "background"  vertical  diffusion 

Large-Scale  Movement  of  Smoke  and  Dust 

■  Use  of  separate  tracer  transport  models 

■  Development  of  simple,  yet  accurate  positive-definite 
tracer  methods  for  use  inside  3-D  models 

Smoke  and  Dust  Removal 

•  Improved  parameterization  of  convective  precipitation 

■  Smoke  and  dust  scavenging  by  precipitation 

■  Smoke  and  dust  effects  on  cloud  microphysics 

•  Smoke  particle  coagulation  in  dense  plumes 

•  Time  evolution  of  chemical  composition  of  smoke 

■  Dry  deposition  at  the  surface 
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NCAR  Global  Model 


Model  Type: 


■  Global  atmospheric  general  circulation  model 
Resolution: 

■  4.5°  latitude  by  7.5°  longitude  grid. 

■  9  layers  in  the  vertical. 

Physical  Processes: 

■  Smoke  transport.  (Spectral  advection) 

■  Smoke  solar  absorption  and  scattering. 

(Fixed  diurnal  mean  Sun) 

■  Smoke  infrared  absorption. 

■  No  smoke  removal. 

NCAR  (1985) 
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Smoke  Scenarios 

Smoke  Mass: 

■  180  Tg  (National  Research  Council,  1985) 

-  60  Tg 

Geographic  Area: 

■  Initially  over  portions  of  NATO  and  Warsaw  Pact 
Height: 

■  Initially  0-7  km. 

Duration: 

■  Smoke  injection  at  a  constant  rate  over  2  days 
No  removal. 

Season: 

■  Northern  Hemisphere  summer  (July) 

■  Northern  Hemisphere  winter  (January) 


NCAR  ( 1985) 


Radiative  Transfer  Enhancements 

Smoke  Scattering  of  Solar  Radiation 
(With  V.  Ramaswamy) 

■  (5  — Eddington  for  Visible  and  Near-IR 

■  Fixed  smoke  optical  properties 

■  Visible  specific  extinction:  5.5  m2g_1 

■  Visible  single  scattering  albedo: 
Low=0.45,  Medium=0.60,  High  =  0.75 

Smoke  Thermal  Infrared  Opac  ity 
(Model  enhancements  by  Curt  Covey) 

■  Gray  absorber  approximation 

■  Specific  absorption  (nominal):  0.5  m2g 
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Conclusions  From  Scattering  Sensitivity  Tests 


For  Fixed  July  Diurnal  Mean  Sun: 

■  Plausible  range  of  scattering  has  little 
effect  on  rate  of  smoke  lofting  and  transport. 

■  For  fixed  specific  absorption, 

going  from  non  -scattering  to  mid-range 
scattering  decreases  land  surface 
temperature  by  about  2°C. 

■  For  fixed  specific  extinction, 

the  plausible  range  of  scattering  changes 
land  surface  temperature  by  ±3°C. 


NCAR  (  1985) 


Surface  Temperature  Over  Land 


Conclusions  From  IR  Sensitivity  Tests 


July  180  Tg  Cases: 

With  IR  Smoke  Opacity  Vs.  Without 

■  Nominal  smoke  IR  opacity  mitigates  land 
surface  temperature  drop  by  roughly  25%. 

■  VERY  large  smoke  IR  opacity  mitigates  land 
surface  temperature  drop  by  roughly  50%. 


NCAR  (1985) 
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Infrared  Effects  of  Saoke : 
Parameterization*  for  Dynamical  Models 


Cnrt  Covey 

RSMAS/MPO.  D.  of  Miami 

J.  T.  Kiehl,  Starley  L.  Thompson 
NCAR 


DNA  Global  Effects  Program  Technical  Meeting 
NASA  Ames  Research  Center,  2S-27  February  1986 


Conclnsions  for  NAS  "baseline"  smoke  amonnts  and  optical  properties : 


1.  In  GCM  calculations,  IR  greenhouse  effect  moderates  "nuclear  winter" 
but  comes  nowhere  near  reversing  it. 


2.  "Gray  smoke"  approximation  provides  a  slight  overestimate  of  the  mag¬ 
nitude  of  the  IR  effect,  compared  with  frequency-dependent  smoke  IR 
absorption  and  emission. 


3.  IR  effects  may  be  more  important  for  mesoscale  modeling  of  initial 
smoke  injection.  The  parameterization  discussed  here  could  be  used 
for  such  models. 


e  =  (oT4)-1  /  (1  -  T J  Bu  du 

a  =  (4oT* )_1  /  (I  -  TJ  (dBw/dT)  du 

where  T  =  temperature 
(i)  =  frequency 

=  blackbody  monochromatic  irradiance 
Tw  =  monochromatic  transaissivity 

=  TH20TC»2T03Ta 

Tx  <=  aonochroaat i c  tranaaisaivity  of  species  z 
("a"  denotes  smoke  and  dnst  aerosol) 


e  and  a  can  be  rewritten  in  the  form 
e  =  eH20  +  *A-H20  +  803Ta(l)  +  ®002Ta ( ’ * * 

where  T  ^  denotes  average  transaissivity  of  x  in  a  frequency  band, 

eH20*  e CO 2 '  e03  *re  eaissivitiea  of  single  species  in  the 
absence  of  other  species,  and 

eA-H20  =  («T4)  /  (1  -  Ta)  TH20  Bw  du 
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FOUR-BAND  PARAMETERIZATION  FOR  e 


A-H20 


4 

eA-H20  =  J  <1  "  Ta(i))TH20(l)p(l).  where 
i=l 


the  four  bands  are  4-8  pm  (  i =1 ) 

8-12  pm  ( i=2) 
12-16  pm  ( i=3 ) 
16-500  pa  ( 1*4) 


T.(i)  =  transmissivity  of  smoke  and  dust  aerosol 
I-(c8'i)ng  +  cd(i'ud)] 


=  exp 


TH20  *  transmissivity  of  water  vapor 
=  Ai  ®xp  (-Bi°e20^ 

p(i)  *  pi«nc^  weighting  function 

L<i>  -  n'i), 


=  Cl' 


C2' (300  -  T) 


u8  =  smoke  column  density  [g  cm*] 
ud  =  dust  column  density  [g  cm  *] 
c  ^  =  smoke  IR  absorption  coefficient  [cm*  g"1] 
cd^  «  dust  IR  absorption  coefficient  [cm*  g_1) 

**1120  =  pressure-weighted  water  vapor  column  density  [g  cm-*] 
T  =  temperature  [I] 


Note : 

Parameters  for  smoke  and  dust  IR  absorption  are  obtained  from  Ramaswamy  and 
Kiehl  (1985).  Parameters  for  water  vapor  IR  absorption  and  Planck  weights 
are  chosen  to  fit  the  narrow  band  computations  of  Rodgers  and  Valshaw 


s', 


~Tempexa} are  dependence. 

(as=  2  g  m'2 ;  4a=o;  aWJo-  ioh3^2) 


T  =  iso  k 


3ho  k 


0.HH51 


0.647? 


'A-U2j0 


inlctyuxL 

aim 

0.6150 


I  So  K 


3ho  k 


^  ~  bcw\ti 

0,503i 
0. 157 1 


0.1%o3 

QMo<Z 


Concj -  basui  (Xpfuioic  foilGu/i.  T- 

oj  Rojn  1  kjjJkL  fo  15%  Caka. 

Oja*  UMlrejnjL  ~f  fCunc^ji, 

BUT  ;  'TIaz,  I flle^iaJL,  dorrriudia,ti<>K.  (does  not  5qJzjl. 

i  nlo  (XcAjuuJbJ:  5kj^  T -  djip^djjAcjL  cj  H9.0 
oJo$o£p  to**  (-for  Cons  is+tncxj  Oji+k  CCM) . 

T- (dep^dd&tGjL  Qjgvaza.  OvhjQ  frorr\  Pla^ruiM 


aged  S 


NCAR  Global  Modeling  1985-1986 


Starley  L.  Thompson 
and 


Filippo  Giorgi 
V.  Ramaswamy 
Curt  Covey 


NCAR  Global  Modeling:  1985-86 


Radiative  Transfer  Enhancements  for  Aerosols 

■  Solar  scattering  and  IR  absorption 

■  Interactive  optical  properties 

Aerosol  Transport  and  Microphysics 

■  Multiple  tracers 

■  Scavenging  parameterizations 


Injection  Scenario  Sensitivity 


NCAR  Global  Model 


Model  Type: 

■  Global  atmospheric  general  circulation  model. 
Resolution: 

■  4.5°  latitude  by  7.5°  longitude  grid. 

■  9  layers  in  the  vertical. 

Physical  Processes: 

■  Multiple  aerosol  transport.  (Spectral  advection) 

■  Aerosol  solar  absorption  and  scattering. 

(Fixed  diurnal  mean  Sun) 

■  Aerosol  infrared  absorption. 

■  Aerosol  scavenging. 


Transport  and  Microphysics  Enhancements 


Radiative  Transfer  Enhancements 


Aerosol  Scattering  of  Solar  Radiation 
(With  V.  Ramaswamy) 

■  <5-Eddington  for  Visible  and  Near-IR 

■  Weighted  5- Eddington  parameters  for  multiple 
aerosol  scatterers 

■  Fixed  smoke  and  dust  optical  properties 

■  Smoke  optical  properties  parameterized  by 
particle  radii.  (Fits  to  Mie  calculations) 

Aerosol  Thermal  Infrared  Opacity 
(Model  enhancements  by  Curt  Covey) 

■  Gray  absorber  approximation 

■  4-band  approximation: 

4-8,  8-12,  12-16,  16-500  /xm 


Aerosol  Optical  Properties:  Solar 


Smoke 

■  "Light" 

■  "Medium" 

■  "Dark" 


m=1.5-0.05i  (cj0 
m=1.5-0.10i  (cj0 
m— 1 .5  -  0.45i  (cjc 


0.75) 

0.60) 

0.45) 


■  Fixed  rg=0.2/xm,  OR  Interactive  radii 
Dust 

■  "Glassy"  melt- derived  basaltic 
(A=0.51/zm,  m=1.53-  O.OOli) 
(A=1.55/zmf  m=1.53-0.002i) 

■  Fixed  rg=0.10/2m  ("Baseline") 

Fixed  rg=0.25/zm  ("Basaltic") 


Aerosol  Optical  Properties:  Infrared 


Non  -  scattering 

Specified  Specific  Absorption 

Smoke 

■  "Grey",  specific  absorption  =  0.5  m2t 

■  Four  broad  bands 

Dust 

■  "Grey",  specific  absorption  =  0.1  m2g 


Dust  Scenarios 


Injected  Dust  Mass: 

■  40  Tg  (National  Research  Council,  1985) 

■  88  Tg  ("Excursion") 

Geographic  Area  of  Injection: 

■  Portions  of  NATO  and  Warsaw  Pact 

Injection  Altitudes: 

■7-15  km.  (Constant  mixing  ratio) 

10  -  24  km.  (Excursion) 

Duration  of  Injection: 

■  Constant  rate  over  1  day. 
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Averaged  Planetary  Albedo 


Smoke  Mass 


Dust- Associated  Suppression  of  Smoke  Lofting 


Fraction  of  Smoke  Above  10  km 
On  Day  15 


Case 


Percentage 


No  Dust 

40  Tg  Basaltic,  rg=0.25/zm 
40  Tg  Basaltic,  rg=0.10/2m 


40 
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Conclusions  From  Scattering  Sensitivity  Tests 


For  Fixed  July  Diurnal  Mean  Sun: 


Plausible  range  of  scattering  has  little 
effect  on  rate  of  smoke  lofting  and  transport. 


For  fixed  specific  absorption, 
going  from  non -scattering  to  mid-range 
scattering  decreases  land  surface 
temperature  by  about  2°C. 


The  range  of  scattering  from  plausible 
variations  in  smoke  composition  changes 
land  surface  temperature  by  ±3°C. 
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Conclusions  From  IR  Sensitivity  Tests 


July  180  Tg  Cases: 

With  Aerosol  IR  Opacity  Vs.  Without 

■  Nominal  aerosol  IR  opacity  mitigates  land 
surface  temperture  drop  by  roughly  25%. 

■  "Grey"  absorber  approximation  for  smoke  is 
reasonably  accurate  if  8-12  /zm  "window" 
specific  absorption  is  used  for  "grey"  value. 


Simplified  History  of  Temperature  Effects 


■  NCAR  Models  With  180  Tg  Smoke  Injection 

■  Mid- latitude  NH  Land  Surface  Temperature 

Model  Temperature  Change 


1-D,  Annual  -35°C 

3-D,  July  ... 

Zonal  Smoke,  No  Removal  -18°C 

Transport,  No  Removal  -22°C 

Transport,  No  Removal,  Aerosol  IR  -15°C 

Transport,  Removal,  Aerosol  IR*  -12°C 


t  Days  5-10 
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Title:  Atmospheric  Effects  of  Nuclear  War  Aerosols  in  GCM  Simulations:  Modeling  of 
Smoke  Microphysical  Processes 

By:  Starley  Thompson,  Filippo  Giorgi  (National  Center  for  Atmospheric  Research) 

A  comprehensive  aerosol  scavenging  package  that  includes  wet  removal,  coagulation,  sedimen¬ 
tation  and  dry  deposition  has  been  incorporated  into  a  global  atmospheric  general  circulation  model 
(GCM)  and  has  been  used  in  simulations  of  nuclear  war  smoke  and  dust.  The  parameterization  uses 
a  new  computationally  efficient  technique  based  on  predicting  the  parameters  of  an  assumed  par¬ 
ticle  size  distribution  rather  than  dividing  the  particles  in  many  finite  size  bins.  Initial  simulations 
assume  a  lognormal  size  distribution  for  smoke  particles  and  a  fixed  particle  radius  dispersion  about 
the  mode  size.  Thus,  only  two  dependent  quantities  are  needed  to  completely  specify  the  smoke 
size  distribution,  e.g.,  smoke  total  mass  and  mode  radius.  Dust  is  treated  simply  as  a  distribution 
having  fixed  mode  radius  and  size  dispersion.  Simulations  have  been  performed  to  determine  the 
sensitivity  of  the  smoke  removal  rate  to  initial  smoke  injection  radius  and  size  dispersion.  The 
evolution  of  smoke  radii  as  a  function  of  time  and  space  are  of  particular  interest  since  the  optical 
properties  of  smoke  depend  on  particle  size.  Results  of  simulations  having  full  feedback  between 
smoke  radii  and  solar  optical  properties  will  be  presented. 


ATMOSPHERIC  EFFECTS  OF  NUCLEAR-WAR  AEROSOLS 

IN  GCM  SIMULATIONS: 

MODELING  OF  SMOKE  MICROPHYSICAL  PROCESSES 


GENERAL  REQUIREMENTS 

•  Comprehensive  aerosol  microphysics. 

A  model  in  which  the  properties  of  the  particle  population 
are  prescribed  cannot  describe  the  effect  of  important  mi¬ 
crophysical  processes. 

•  Computational  efficiency. 

A  detailed  bin  model  is  too  expensive  for  use  in  GCMs. 


•  Adopt  an  intermediate  approach. 


MODEL-AEROSOL  CHARACTERISTICS 


•  The  model-aerosol  consists  of  modes,  each  characterized 
by 

-  Homogeneous  chemical  composition 

-  Log-normal  size  distribution 


-(In r  -  lnrg)2' 


dlnr 


The  defining  parameters  of  the  log-normal  form  are 
N  =  Total  Particle  Number  Concentration 
r9  =  Geometric  Mean  Radius  (by  number) 

<x  =  Geometric  Standard  Deviation 


•  Why  the  log-normal  form  ? 

-  Representative  for  ambient  aerosols 

-  Convenient  analytical  and  numerical  properties 


•  Assume  that  the  evolution  of  the  particle  size  distribution 
can  be  described  (approximately)  by  the  evolution  of  the 
free  defining  parameters  of  the  corresponding  log-normal 
modes. 


Solve  prognostic  equations  for  an  adequate  number  of  mo¬ 
ments  of  the  log-normal  modes.  The  moment  of  order  j  is 
defined  by 

J  r3n(r)dr 

Make  use  of  the  analytical  relation 


to  obtain  N.  ^  ,  (T^ . 


For  nuclear  smoke  assume  a  fixed  <7g  =  0.3;  free  pa¬ 
rameters  are  N  and  rg  .  Solve  an  equation  for  the  two 
moments 


N  =  J  n(r)dr 
V  =  \wf  r3n(r)di 


Total  number  concentration 


Total  volume 


For  nuclear  dust  assume  a  fixed  crg  =  0.7  and  a  fixed 
rg  =  0.25  fim.  There  is  only  one  free  parameter,  thus 
solve  only  an  equation  for  V . 


•  This  method  requires  integrations  over  the  size  distribution 
of  size-dependent  terms. 

-  The  integration  is  performed  analytically  if  the  term 
can  be  expressed  as  a  polynomial  in  the  particle  radius 
r. 

-  Otherwise  the  integration  is  performed  via  a  computa¬ 
tionally  efficient  technique  that  makes  use  of  a  moving 
”  Hermite-quadrature  grid”  centered  around  rq  . 


•  The  following  processes  are  included  in  the  model: 


-  Sources  (different  scenarios) 

-  3D  Transport  (CCM  — ►  large  scale  advection,  sub¬ 
grid  scale  eddy  diffusion) 

-  Coagulation 

-  Sedimentation 

-  Dry  Deposition 

-  Wet  Removal 

•  The  prognostic  equation  for  the  moments  of  the  particle 
modes  has  the  form 


dN,V 

dt 


+ 


dN,V' 

+ 

aourc 

dN,V' 

+ 

transp 

dN,  V 

dt 

dt 

dt 

i  } 

+ 

atdim 

dN,V' 

dt 

+ 

drydep 

dN,  V 

dt 

+ 


coag 


wetrem 
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COAGULATION 


At  the  present  stage  this  term  includes  only  intramodal 
Brownian  Coagulation 


dN 

dt 


j  coag 


~ 2KnN 2 


.where 


dV 

dt 


=  0 


coag 


f  CXD(  (lnr  -  lnrg)2  \  /  (Ins-  In  r„)2 

/  2 TO*  eXP  l  2 a]  yeXPV  2 a] 


and  K(r ,  5)  is  the  Brownian  Coagulation  kernel 


SEDIMENTATION 


•  This  term  is  included  as  the  vertical  divergence  of  the  grav¬ 
itational  settling  flux. 

diV ,  V  _  dixy 

.  dt  J  stdim  dz 

where  I^Y  are  the  integrated  particle  gravitational  settling 
fluxes 

In  =  J  vg(r)n(r)dr 

Jv  =  “Tr  J  r3vg(r)n(r)dr 

and  tj '  is  the  particle  gravitational  settling  velocity. 

•  The  vertical  derivatives  are  approximated  via  a  mass  con¬ 


serving  finite  difference  scheme. 


DRY  DEPOSITION 


The  dry  deposition  term  is  calculated  at  the  bottom  model 
level  as 

\dN.V\ 

ex  —Dfijy 


dt 


J  drydep 


where  J0Wvare  the  integrated  particle  dry  deposition  fluxes 


and  Vj  is  the  particle  deposition  velocity. 
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The  particle  deposition  velocity  is  parameterized  in  terms 
of: 

-  Surface  momentum  drag  coefficient 

-  Wind  speed,  air  density,  and  temperature  at  the  bot¬ 
tom  CCM  level 

-  Particle  radius  and  density 

-  Characteristics  of  the  surface  roughness  elements 

Deposition  processes  included  in  the  parameterization  are: 

-  Turbulent  transport  through  the  Boundary  Layer 

-  Brownian  Diffusion,  Interception,  and  Impaction  on 
the  surface  roughness  elements 

Surface  types  considered  are 

-  Land  (treated  as  a  grassy  surface) 

-  Ocean 

-  Snow  and  Sea  Ice 


WET  REMOVAL 


•  The  model  includes  only  in-cloud  scavenging. 

•  The  assumption  is  made  that  the  aerosol  is  hygroscopic 
and  that  inside  clouds  at  the  onset  of  precipition  it  resides 
entirely  in  the  liquid  phase.  The  local  particle  removal  rate 
is  thus  equal  to  the  liquid  water  removal  rate  A  (s'1),  i.e. 

=  -\{N,V] 

wetrem 

(Note  that  this  term  does  not  affect  the  value  of  r3 ). 


'dN,y 

dt 


A  is  calculated  in  terms  of  the  local  CCM-produced  conden¬ 
sation  rate  Q  as  follows: 

-  Since  precipitation  is  a  sub-grid  process  the  following 
quantities  are  first  defined: 

F  =  fraction  of  the  grid  box  where  precipitation  oc¬ 
curs 

/ 3  =  in-cloud  liquid  water  removal  rate 
L  =  in-cloud  liquid  water  content 
Tb  =  duration  of  the  precipitation  event  within  the 
time  step  At. 

Q  -  FLt% 


-  At  each  grid  point  and  time  step  A  is  given  by 


X  =  ~  fe-^-ll 
At  1  J 


-  F.  P,  L,  Tb  are  expressed  in  terms  of  Q  and/or  given 


parameters.  Different  parameterizations  are  used  for 
stable  and  unstable  condensation. 


•  Validation  of  wet  removal  parameterization. 


-  CCM-derived  TROPOSPHERIC  RESIDENCE  TIMES 
(t)  of  transparent  aerosols  with  different  source  dis¬ 
tributions: 

Surface  injection  — ►  T  =  2-3  days 
Lower  tropospheric  injection  (z<7  km)  — * 

T  =  5-7  days 

Middle  tropospheric  injection  (4<z<10  km)  — * 

T  =  10-15  days 

Stratospheric  injection  — ►  r  =  25-35  days 

-  Generally  good  agreement  is  found  with  available  es¬ 
timates  from  atmospheric  radionuclide  observations. 


EFFECTS  OF  MICROPHYSICS  IN  INTERACTIVE 
NUCLEAR  WINTER  SIMULATIONS 


Sedimentation  and  Dry  Deposition 

-  No  major  effects  in  30-day  simulations. 

-  Possibly  more  important  on  a  longer  time  scale  or  for 
hydrophobic  aerosols. 

Wet  Removal 

-  Very  important  because  it  determines  smoke  atmo¬ 
spheric  residence  times. 

-  The  smoke  residence  time  increases  from  a  value  of 
the  order  of  a  few  days  for  the  case  of  transparent 
aerosols  to  a  value  of  the  order  of  100  days  when  the 
radiative  effects  of  smoke  have  altered  the  thermal  and 
dynamical  structure  of  the  atmosphere. 

Coagulation 

-  Of  relevance  because  it  affects  the  value  of  the  smoke 
geometric  mean  radius  and.  in  turn,  the  smoke  bulk 
optical  properties. 
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Radius 


y  Averaged  Smoke  Radius 
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OUTLINE 


Photochemical  and  Climatic  Effects 
of  Trace  Gases  Injections 

•  Previous  studies 

o  Estimated  inputs  compared  with  current 

global  burdens  and  fluxes 

©  Effects 

Direct  Climate 

Indirect  O3 ,  tropospheric  OH 

climate,  chronic 
biological  impact, 
military  operations 

o  Interactions — relations  to  other  areas 


Previous  Studies  (Pre-1980) 
Effect  of  NO  on  0 

x  3 

Effects  of  oxides  of  nitrogen  on 
stratospheric  O3 

NO  +  O3  4*  NO2  +  O2 
NO2  +  0  >  NO  +  O2 


Tropospheric  ozone 


R00  +  NO  -►  N02  +  RO 
NO2  +  hv  NO  +  0 

0  +  O2  ^  O3 


higher  NOx  =>  higher  O3  production 


Net  effect  of  N0X  injections  depend  on 
amount  and  altitude  distributions  and  R00 


Trace  Gas  Inputs 


Trace  Gas  Inputs 


Direct  Effect 


NO2 


Solar  Heating 


N2O 

H2O 

Hydrocarbon,  CH4  f 


Greenhouse  Effect 


HNO3 


Dust  and  Smoke 


Obscuration 
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Indirect  Effect 


N02  +  HNO3 

CH4  +  H2O 

CO,  CH4  ■>  suppress  OH 
Dust  and  smoke  heterogeneous  chemistry 


Climate 


Addition  of  greenhouse  gases 


Redistribution  of  O3 


Photochemical 


Stratosphere 

O3 


Troposphere 
OH,  hydrocarbons,  O3 


Ozone  Chemistry 


Production 

Photolysis  of  O2 
O2  +  hv  -►  0  +  0 

0  +  02  ^  03 

Removal 

Oxygen 

H0X 

N0X 

ciox 


Smog  Chemistry 
R00  +  NO  +  NO2  +  RO 
NO2  +  hv  +  0  +  NO 
0  +  O2  O3 

0  +  O3 

OH  +  0;  0  +  HO2 

OH  +  O3 ;  O3  +  HO2 

NO  +  03  +  NO2  +  02 
NO2  +  0  +  NO  +  O2 

Cl  +  O3  +  CIO  +  O2 
CIO  +  0  *  Cl  +  02 


Things  that  affect  ozone  (ppm) 


CH4, 


N2O 


CFCs 


H2O  (ppm) 

OH,  HO2  (ppb) 


(ppb) 

NO,  N02,  HNO3  (ppb) 


(ppb) 

Cl,  CIO,  HC1  (ppb) 


Tropospheric  OH 


Production 


03  +  hv  >  O^D)  +  02 
O^D)  +  H20  >  20H 


Removal 


OH  +  CH4 

OH  +  CO- 


14K 


Cleansing  agent  in  the  troposphere 
controls  concentration  of  gases  such 
as  (CH3 )2S ,  H2S,  hydrocarbons,  CH2CI 

The  lifetimes  of  some  of  these  gases 
is  ~  days 

If  production  of  these  gases  is  main 
tained  while  OH  is  suppressed,  the 
concentrations  will  increase  signifi 
cantly  within  days. 


No  sunlight,  no  OH 

Build  up  of  (CH3)2S,  H2S,  hydrocarbons 

OH  production  resumes  after  return  of  sun¬ 
light;  however,  CO,  CH4  concentrations  are 
much  higher.  OH  concentration  could  be 
much  smaller.  Recovery  time  will  take 
years. 

Question :  Could  enhanced  penetration  of  UV  radi¬ 
ation  due  to  removal  of  stratospheric 
O3  and  higher  tropospheric  concentra¬ 
tion  of  O3  from  smog  reactions  compen¬ 
sate  for  the  increased  CO  and  CH4? 
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Ozone  and  other  processes 

Heat  budget  radiation  circulation 

Optical  properties  of  atmosphere 
photolysis  rates 

Chemistry  O^D)  +  H2O  H0X 

O^D)  +  N20  NOY 

Shielding  of  UV  radiation 

Smog  chemistry 


Effect  on  plants 


Effects  of  Transport  on  Ozone 


Zonal-mean  Picture 


Chemistry 


Composition 


»J 


AER  MODELS 


1-D  photochemical-diffusion  model 


1-D  diurnal  model 


1-D  radiative  convective  photochemical  model 


2-D  zonal-mean  model 

o  diabatic  circulation  formulation 

•  classical  Eulerian  formulation 


2-D  climate  model 


3-D  semispectral,  global  primitive  equation 
model 
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Conclusion 


Atmospheric  radiation,  dynamics,  and 
chemistry  are  highly  coupled 


Identify  areas  of  interest  using 
simple  models 


Solutions  of  problems  through  coupled 
models 
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ABSTRACT 

In  a  nuclear  winter  scenario,  smoke  and  soot  fall  out  of 
the  atmosphere  and  lower  the  albedo  of  the  snow  and  ice.  We  use 
an  energy  balance  climate  model  to  test  the  effect  this  would 
have  on  moderating  the  nuclear  winter  cooling.  We  find  that  the 
effect  is  small,  since  the  maximum  effect  on  surface  albedo 
occurs  when  the  atmosphere  still  contains  a  large  amount  of 
soot,  thereby  reducing  the  importance  of  surface  albedo. 

The  climate  model  illustrates  that  snow  and  ice  feedbacks 
can  prolong  the  surface  temperature  effects  of  nuclear  winter 
for  several  years. 


TEXT 

Following  a  nuclear  war,  burning  cities  would  put  smoke  and 
soot  into  the  atmosphere.  The  amount  of  smoke  and  its  optical 
properties  depend  on  many  factors,  as  discussed  in  other  papers 
in  this  volume.  When  the  smoke  is  removed  from  the  atmosphere, 
it  will  lower  the  albedo  of  the  snow  and  ice  upon  which  it  falls 
(Warren  and  Wiscombe,  1985) .  In  this  paper  we  use  several 
hypothesized  nuclear  winter  scenarios  and  the  albedo  effects 
described  by  Warren  and  Wiscombe  to  calculate,  using  an  energy 
balance  climate  model,  the  effects  of  dirty  snow  on  moderating 
the  cooling.  The  figure  captions  discuss  the  results  in 
detail.  In  this  text,  the  experiments  we  performed  are  briefly 
described.  A  longer  paper  describing  the  research  in  detail  is 
in  preparation. 

The  climate  model  of  Robock  (1983,  1984)  was  used  for  these 
experiments.  We  modified  the  results  of  Warren  and  Wiscombe 
(1985)  to  give  clean  snow  and  ice  albedos  to  correspond  to  the 
values  used  in  the  model's  surface  albedo  parameterization  (Fig. 
1)  • 

Baseline  nuclear  winter  scenarios  from  Turco,  _ aJL^. 

(1983)  (TTAPS) ,  National  Academy  of  Sciences  (1985)  (NAS),  and 
Malone,  et  al.  .  (1986)  were  used  to  force  the  climate  model.  The 
atmospheric  smoke  mass  loading  in  each  of  these  cases  is  225  Tg 


(TTAPS) ,  180  Tg  (NAS)  and  170  Tg  (Malone,  et  al  .) .  A  radiative 
equilibrium  smoke  cloud  model  calculated  the  downward  shortwave 
and  longwave  radiation  from  the  cloud  as  a  function  of  optical 
depth  at  each  latitude  (Crutzen,  et  al . .  1984) .  For  scattering 
calculations,  the  assymetry  factor  was  held  constant  at  0.7 
while  the  single  scattering  albedo  is  varied  to  produce  the  same 
solar  insolation  as  published  in  the  NAS  and  TTAPS  papers.  We 
assumed  the  infrared  optical  depth  was  one  tenth  the  optical 
depth  of  the  visible  radiation  and  ignored  scattering  of 
infrared  radiation. 

The  soot  component  of  the  smoke  was  used  to  calculate  the 
snow  and  ice  albedos  (Warren  and  Wiscombe,  1985) .  From  the 
information  in  the  above  papers,  we  calculated  the  fallout  rate 
of  soot  onto  snow  as  a  function  of  time  and  latitude.  Dirty 
snow  calculations  for  the  Malone  et  al .  case  were  not  carried 
out  because  we  could  not  obtain  information  about  the  soot 
fallout.  In  order  to  calculate  mass  fraction  of  soot  in  snow, 
we  assumed  a  precipitation  rate  of  1  cm  liquid  water  per  month. 
This  is  a  very  low  rate  for  the  current  climate,  but  was  chosen 
as  a  "worst  case."  Some  simulations  suggest  reductions  of 
precipitation  as  a  result  of  nuclear  winter.  GCM  calculations 
of  snowfall  and  soot  mass  fraction  were  kindly  provided  by  Steve 
Ghan  (Ghan,  et  al . .  1985) .  Our  values  give  the  same  or  higher 
soot  mass  fractions  than  theirs  and  the  same  or  lower 
precipitation  rates.  Any  smaller  mass  fractions  of  soot  would 
give  less  of  an  albedo  effect.  We  ignore  for  now  the  possible 
reappearance  of  soot  after  melting  on  ice  caps  or  sea  ice,  which 
was  suggested  by  Warren  and  Wiscombe  (1985) . 

We  ran  the  model  for  the  various  nuclear  winter  scenarios 
with  and  without  dirty  snow,  and  calculate  the  differences  in 
the  surface  temperature  response.  The  largest  effects  are  for 
the  TTAPS  scenario,  and  these  results  are  shown  in  detail.  We 
also  demonstrate  the  model  sensitivity  to  snow/albedo  and  sea 
ice/thermal  inertia  feedbacks  separately. 

Dirty  snow  is  found  to  have  small  effect  on  ameliorating 
the  cooling  induced  by  massive  soot  clouds  in  the  atmosphere. 
The  largest  effect  is  for  a  soot  injection  in  spring,  allowing 
the  atmosphere  to  be  clean  enough  in  summer  at  high  latitudes 
for  the  large  insolation  to  interact  with  the  darker  snow. 
Still,  in  this  case,  the  effect  is  only  about  20%  of  the  total 
cooling  at  the  most. 
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FIGURE  CAPTIONS 


1 .  Albedo  of  snow,  sea  ice,  and  old  snow  and  sea  ice  as  a 

function  of  soot  mass  fraction  in  snow.  See  Robock  (1983) 

for  snow  and  ice  albedo  parameter izations  and  Warren  and 
Wiscombe  (1985)  for  soot  effects. 

2.  Response  of  model  to  lowering  albedo  of  snow  and  ice  by 

50%.  Shown  are  surface  air  temperature  differences  from  an 
unperturbed  reference  run  (°C)  for  first  three  years 
following  change  (a)  and  at  equilibrium  (b)  as  functions  of 
latitude  and  time  of  year.  The  largest  effect  is  in  the 

summer  at  high  latitudes. 

Extinction  optical  depth  of  smoke  as  a  function  of  time  for 
t.he  scenarios  considered.  For  TTAPS  and  NAS,  the  smoke  is 
spread  uniformly  throughout  the  Northern  Hemisphere.  For 
the  Malone  et  al .  cases,  the  latitudinal  spread  is  time 
dependent  as  given  in  the  paper.  For  times  after  the  end  of 
the  GCM  run,  the  smoke  is  assumed  to  not  spread  any  more  and 
to  decay  exponentially  in  time  at  the  rate  specified  in  the 
paper.  For  comparison,  an  area  weighted  average  was  taken 
of  these  optical  depths  to  obtain  a  Northern  Heimspheric 
average . 


Albedo  of  new  snow  as  calculated  for  the  TTAPS  and  NAS 
cases.  We  did  not  have  enough  information  about  soot 
fallout  rates  to  do  this  calculation  for  the  Malone  et  al . 
case . 

Surface  temperature  changes  (°C)  for  the  TTAPS  summer  case 
for  the  first  5  years  of  the  response  (differences  from 
reference  run) .  High  latitude  enhancement  of  the  response 
is  evident. 

Demonstration  of  the  cryosphere  feedbacks.  a.)  Difference 
in  response  for  a  case  with  no  snow  and  ice  feedbacks  (fixed 
areas)  and  control  (with  all  feedbacks,  as  shown  in  Fig. 
5) .  b.)  As  in  a.,  but  with  only  snow  area  responding.  The 
summertime  enhancement  is  evident.  There  is  little 
snow/albedo  feedback  in  winter  when  the  sun  is  weak.  c.) 
As  in  a.,  but  with  only  sea  ice  area  responding.  The 
maximum  effect  is  in  the  winter  at  high  latitudes,  because 
the  changing  sea  ice  changes  the  thermal  inertia  affecting 
the  amplitude  of  the  seasonal  cycle.  This  feedback 
dominates  the  cryosphere  feedback. 

As  in  Fig.  5,  but  for  NAS  scenario.  The  response  is  smaller 
and  does  not  last  as  long. 

As  in  Fig.  5,  but  for  Malone  et  al .  case.  The  initial 
response  is  smaller  than  the  NAS  case  (Fig.  6),  but  is 
longer  lasting  because  of  the  longer  soot  residence  time 
(Fig.  3)  . 

For  a  TTAPS  scenario  with  dirty  snow  (soot  case)  starting  in 
winter,  the  surface  temperature  response  (a),  the  surface 
temperature  difference  between  the  soot  case  and  control 
case  (b) ,  and  the  albedo  difference  (c) .  The  cooling  is 
about  1®C  less  with  dirty  snow  than  without  it  at  high 
latitudes  in  the  summer  (b) .  The  largest  albedo  changes  due 
to  dirty  snow  are  in  the  winter,  but  the  maximum  effect  is 
in  the  summer,  due  to  both  a  cleaner  atmosphere  and  greater 
insolation  then. 

As  in  Fig.  9,  but  for  spring  for  3  years.  This  case  shows 
the  largest  dirty  snow  effect,  of  more  than  2°C  the  first 
summer  and  1°  C  the  second  summer . 

As  in  Fig.  9,  but  for  summer  for  3  years.  The  actual 
cooling  is  greatest  for  this  case,  and  the  dirty  snow 
effect  is  smaller.  This  is  because  by  the  time  the 
atmosphere  becomes  clean  enough  for  surface  albedo  to 
matter,  it  is  already  fall,  and  insolation  has  decreased. 
The  fall  case  has  the  smallest  actual  response  (Fig.  15) 
and  dirty  snow  effect,  not  shown. 


A  five  year  run  of  the  TTAPS  case  with  dirty  snow  for 
winter . 


A  five  year  run  of  the  TTAPS  case  with  dirty  snow  fo 
spring . 

A  five  year  run  of  the  TTAPS  case  with  dirty  snow  fo 
summer . 


A  five  year  run  of  the  TTAPS  case  with  dirty  snow  fo 
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AN  EXPERIMENTAL  PROTOCOL  FOR  AIRCRAFT  OBSERVATIONS  OF  THE  AEROSOLS  AND  GASES 


FROM  LARGE  FIRES 

Lawrence  F.  Radke 
Cloud  and  Aerosol  Research  Group 
University  of  Washington 
Seattle,  WA  98195  U.S.A. 

The  robustness  of  the  nuclear  winter  scenario  depends  critically  on  a 
number  of  uncertain  factors,  some  of  which  can  be  experimentally  examined  by 
studying  large  prescribed  fires.  We  describe  here  the  experimental  protocol 
used  to  study  the  following  Issues: 

1)  "Prompt  removal"  of  smoke  by  fire  capping  cumulus  clouds 

2)  Unusual  microphysical  characteristics  of  both  fire  capping  cumulus  and 
other  clouds  later  Ingesting  the  smoke. 

3)  Variations  in  particle  and  gaseous  emissions  with  fire  size  and  fuel  mix 

4)  Radiative  properties  of  the  smokes 

5)  Time  evolution  of  dense  smokes  with  emphasis  on  aerosol  coagulation, 
and  changes  in  morphology  and  radiitive  properties. 


Abstract  for  the  25-27  February  1986  Global  Effects  Progran  Technical  Meeting. 

Philip  J.  Riggan 
USDA  Forest  Service 

Pacific  Southwest  Forest  and  Range  Experiment  Station 

4955  Canyon  Crest  Drive 

Riverside,  CA  92507 

FTS  796  6516/  (714)  351  6516 

title:  The  1986  field  experiment  at  the  San  Dimas  Experimental  Forest  (Lodi 
Canyon) ,  Los  Angeles  County. 

An  overview  will  be  presented  of  the  current  plans  for  the  420  hectare,  high- 
intensity  prescribed  fire  in  Los  Angeles  County.  Plans  for  conducting  the 
fire,  research  coordination,  expected  fire  behavior,  measurements  of  fire 
physics,  and  organization  of  fire  management  will  be  discussed. 


•  EMISSIONS  FROM  OPEN  FIRES  ARE  KNOWN  TO  BE  HARMFUL  TO 

HUMAN  HEALTH. 

•  EMISSIONS  FROM  PRESCRIBED  FIRES  CAN  BE  REDUCED. 

•  SOURCE  STRENGTH  DATA  ARE  NEEDED  TO  MANAGE  SMOKE. 

•  FINGERPRINTING  METHODS  ARE  NEEDED  TO  IMPROVE  THE 

ASSESSMENT  OF  IMPACTS  AT  RECEPTOR  SITES. 


WHAT  IS  RESEARCH  DOING? 


□  Techniques  of  burning  to  minimize  the  production  of 

AIR  POLLUTANTS  (SOURCE  MODELS) 

□  Update  national  and  regional  emission  inventories 

(Sandberg  and  Peterson  1985,  Chi  et  al.  1979) 

O  Tradeoffs  between  smoke  from  wild  fires  and 
prescribed  fires  (Cooper  1976) 

□  An  assessment  of  emissions  production  as  a 

FUNCTION  OF  INCREASED  WOOD  UTILIZATION 

(Sandberg  and  Ward  1982) 

a  A  model(s)  to  provide  data  ON  SOURCE  STRENGTH  for 
SMOKE  DISPERSION  CALCULATIONS  (PlEROVICH  1983) 


□  Source  characteristics  for  assessing  impact  of 

BURNING  OF  FOREST  MATERIALS  ON  REGIONAL  HAZE 
(Core  ET  AL.  1983,  1984,  1985) 


a  Database  of  emissions  data  and  research. 
—Need  to  expand  on  effort  underwayw 


••(1986)  WHAT  ARE  WE  DOING  NOW? 


□  Piled  slash  (Westside  of  Cascades)— 6  samples 

□  Broadcast  burns  of  hardwood  slash— 6  samples 

□  Broadcast  burns  of  long-needled  conifer  slash  east  of 

Cascades— 4  to  6  samples. 

□  Lodi  Experiment— 4  to  6  samples  of  chaparral? 

••(  1 987-  )  WHAT  WILL  WE  BE  DOING? 

□  Slash  VS.  Grass 

O  Herbicide  treated  units 

PRODUCTS  OF  DECOMPOSITION 

Effectiveness  of  herbicide  treatments. 

□  Partially  cut  units  (Line  fires  in  conifer  types). 
p  Piled  slash  VS,  Broadcast  VS,  Partially  cut 

□  Brush  species 

Bitterbrush 
Sagebrush 
Ceanothus 
Chaparral 

d  Combustion  and  fuel  chemistry  affecting  production 
of  trace  materials.  (Combustion  hood  research.) 

*  j v 
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A.  Surface  Sampling 

1.  vertical  array— LINE  FIRES 

2.  Horizontal  array— AREA  FIRES 

3.  Traversing  in  horizontal  plane— PILED  MATERIAL  OR 

POOL  FIRES 

B.  Airborne  Sampling  systems 


FMISSION  MEASUREMENTS 


A.  Emission  factors 


1 .  CARB0N-MA55  BALANCE 


2.  Emission  flux/fuel  consumption 


B.  Source  strength 


1.  Emission  flux 


2.  Fuel  consumption  x  EF 


5L 
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SAMPLING  MATRIX 


I.  Particulate  matter 

A.  47-MM  GLASS  FIBER 

B.  47-mm  glass  fiber  backed  with  polyurethane  foam 

C.  37-MM  QUARTZ 

D.  37-mm  teflon 

E.  20-  X  25-CM  GLASS  FIBER  (Hl-VOL) 

II.  Gases 

A.  Grab  samples 

1 .  Bag  samples  for  combustion  gas  analysis. 

2.  Adsorption  tubes. 

3.  Canister  samples  for  trace  gases. 

B.  Real  time  gas  analysis 

1 .  CO  (  I  TO  40  PERCENT  OF  CARBON) 

2.  C02  (50  TO  95+  PERCENT  OF  CARBON) 

III.  OTHER  SENSORS 

A.  Temperature  (Type  K  T/C,  0  #C  compensated) 

B.  Combustion  gas  velocity  (Yane  anemometers,  Kurz  mass- 

flow  meters) 

IV.  Data  logger  (Fluke  2280A) 

Y.  Gas  sampling  is  continuous 

A.  5  TEFLON  TUBES  100  M  X  3/16"  ID  (2  LITERS/MIN) 

B.  3  Nylon  tubes  100  m  x  1/4"  ID  (10  liters/min) 

C.  2  Nylon  tubes  100  m  x  1/4"  ID  (5  liters/min) 


3/4 


IM-AltS  1S1  TECHNICAL  PAPERS  PRESENTED  AT  THE  DEFENSE  NUCLEAR 

AOENCV  GLOBAL  EFFECTS  R.  .  <U)  DOD  NUCLEAR  INFORMATION 
AND  ANAL VS IS  CENTER  SANTA  BARBARA  CA.  .  IS  HAV  SC 
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Particulate  matter 


A.  PM 

1-  Total  Particulate  Matter  EFprgf:]  and  EFpnis] 

2.  POM  profiling  (B[a]P) 

3.  Light  absorption 

B  PM2.5 

1 .  Particulate  matter  less  than  a  mean  cutpoint 

DIAMETER  OF  2.5  JIM  (EFpf-|2.5lFl  AND  EFpn2.5l5l) 

2.  Carbon  analysis 

A.  GRAPHYTIC 

b.  Organic 

c.  Total  carbon 

GASES 


A.  Bag  grab  samples  (CO,  C02,  CH4,  NMHC) 

B.  CANISTERS  (CH3CL,  N02,  CH4,  C02,  COS,  ETC.) 

C.  Adsorption  tubes  (Aldehydes  and  other  heavy  MWHC) 
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IX  Percent  elemental  composition  of  PM2.5  for  specific  trace 
elements  and  ions  by  combustion  phase. 


Table  XII  Percent  elemental,  organic,  and  total  carbon  content 
combustion  phase. 


e  VIII  Correlation  matrix  for  variables  measured  during  1983  and  1984 
tests  (F  and  SI  data  only). 


COne  EFFIC  FLAME  □  COMB  EFFIC  SI  ■  COMB  EFFIC  S2 


1.  WAYCROSS  GE0R6IA  (WARD  ET  AL.  1979) 

a.)  Airborne  (flux  of  PM  through  vertical  planes  at 

SPECIFIED  DISTANCES  DIVIDED  BY  THE  INVENTORIED 
MEASUREMENT  OF  FUEL  CONSUMPTION. 

805  M  Q  =  253.1  G/S  W  =  13.08  KG/S  EFpn  =  19.4  6/KG 

3360  M  Q  =  193.0  G/S  W  =  13.08  KG/S  EFpn  =  14.8  G/KG 

7355 MQ  =  255.8  G/S  W  =  1 3.08  KG/S  EFpn  =  1 9.6  G/KG 

805  M  Q  =  236.2  6/S  W  =  1 4.55  KG/S  EFpn  =  1 6.2  G/KG 

b.)  Surface  towers  (CMB)  EFpm  =  8.1  6/kg  Mass  flux  EFpn  = 

1  1.2  G/KG;  (FLUX  OF  PM  THROUGH  VERTICAL  PLANE  DIVIDED  BY 
THE  SURFACE  MEASUREMENT  OF  FUEL  CONSUMPTION.) 

C.)  CCN  PARTICLE  SIZE  DISTRIBUTION  OVER  DISTANCE.  VOLUME 

FLUX  METHOD.  RELATIVE  HUMIOITY  VARIED  FROM  65%  TO  98  %. 
THE  PARTICLE  SIZE  DISTRIBUTION  CHANGED  AS  THE  PLUME 
MOVED  DOWNWIND.  MASS  FLUX  OF  THE  PM  AS  DETERMINED 
USING  THE  INTEGRATING  NEPHELOMETER  REMAINED  FAIRLY 
CONSTANT  FOR  THE  THREE  TRAVERSES.  THE  PM  FLUX  AS 
DETERMINED  USING  THE  CN  DATA  DROPPED  OFF  RAPIDLY.  THE 
TOTAL  NUMBER  OF  PARTICLES  MOVING  DOWNSTREAM  DECREASED 
BY  A  FACTOR  OF  2.6,  BUT  THERE  WAS  AN  INCREASE  IN  THE 
NUMBER  OF  PARTICLES  BETWEEN  0.3  AND  0.5  JIM.  THE  NUMBER 
OF  PARTICLES  BETWEEN  0.5  AND  10  JIM  DECREASED.  A 
CONSTANT  RESPONSE  FROM  THE  INTEGRATING  NEPHELOMETER 
WHILE  THE  TOTAL  NUMBER  OF  PARTICLES  DECREASED  SUPPORTS 
THE  IDEA  OF  PARTICLE  COAGULATION  IN  TRANSIT.  IN  EFFECT, 
THE  INCREASED  SCATTERING  EFFICIENCY  DUE  TO  THE  PARTICLES 
ACCUMULATING  IN  THE  0.3  TO  0.5-JlM  SIZE  RANGE  IS  THOUGHT 
TO  OFFSET  THE  LOSS  OF  SCATTERING  BY  THE  GREATER  NUMBER 
OF  SMALL  DIAMETER  PARTICLES  OBSERVED  FOR  THE  805-M 
CROSS  SECTION. 


2.  Green  Mountain  Oregon  (Ward  et  al.  1982) 

a. )  Emissions  flux  from  airborne  samples 

(PM:  526, 637.  417,  AND  767  6/5) 

b. )  Rate  of  fuel  consumption  from  ground  measurements. 

(W:  28.1 , 23.4,  20.5,  AND  41 .2  KG/S) 

C.)  EFpfi  VALUES  FROM  AIRBORNE  SAMPLES. 

(EFW:  18.7,  27.2,  20.3,  AND  18.6  G/KG). 
d)  Surface  based  CMB  method  for  determining  EFpn  values. 
(EFpmif]:  12.0  G/KG,  EFpM[S]:  21  G/KG  AND  28.5  6/KG, 
EFpfi  =  0.5*EFpm{f]  +  0.5*EFpfi[sj  =  20.3  G/KG.) 


3.  Joule  Oregon  (Radke  et  al.  1986,  Ward  1986) 

4.  Weyerhaeuser  (Washington)  FS  (Radke  et  al.  1986,  Ward 

and  Sandberg  1 986) 

5.  Weyerhaeuser  (Washington)  U.  of  W.  (Stith  et  al.  1981 ) 
a.)  Particle  volume  distribution  peaked  at  0.3  jim 

B. )  EFpfi  =  0.2  TO  2.0X. 

C. )  THE  BURNS  WITH  THE  HIGHEST  RATES  OF  FUEL  CONSUMPTION 

HAD  THE  HIGHEST  DENSITY. 

D. )  Large  source  of  CCN;  *  IO*o  to  10* 1  CCN  were  produced 

PER  GRAM  OF  WOOD  CONSUMED.  THEY  PRODUCE  UNUSUALLY  HIGH 
CONCENTRATIONS  OF  SMALL  DROPLETS  (<  10pM  DIAMETER). 

e.)  Large  reduction  on  visibility  (in  one  plume  the  visual 

RANGE  WAS  3%  OF  THAT  IN  CLEAN  AIR). 


6.  Australia  (Evans  et  al.  1976) 

a. )  The  particulate  matter  yield  from  a  hot  clearing  burn 

was  reported  to  BE  .44% 

b. )  The  particulate  matter  yield  from  a  fuel  reduction  burn 

WAS  2%  TO  4%. 

c. )  The  photochemically  reactive  species  may  be  5  times 

MORE  CONCENTRATED  IN  THE  HOT  CLEARING  BURN  THAN  IN  THE 
FUEL  REDUCTION  BURN. 

D.)  THE  HOTTER  THE  FIRE  THE  LARGER  THE  PRODUCTION  OF 

NITROGEN  OXIDES.  ALL  OF  THE  AUSTRALIAN  FIRES  WERE  LOW 
INTENSITY  BURNS. 
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A.  TO 


COMPARE  AIRBORNE  MEASUREMENTS  WITH  SURFACE-BASED 
MEASUREMENTS  FOR  HIGHER  INTENSITY  FIRES  (LOTI  CANYON, 
CHABLEAU  TYPE  FIRES). 

B.  TO  COMPARE  MEASUREMENTS  OF  EMISSIONS  FROM  POOL  FIRES 

USING  EQUIPMENT  FOR  MEASURING  EMISSIONS  FROM  FOREST  AND 
RANGE  FIRES. 

C.  TO  MEASURE  EMISSIONS  FROM  FIRES  OF  RUBBLE! ZED  MATERIALS 

FROM  URBAN  ENVIRONMENTS  WITH  EQUIPMENT  USED  FOR 
MEASURING  EMISSIONS  FROM  FOREST  AND  RANGE  FIRES. 

D.  TO  ASSESS  FACTORS  CONTRIBUTING  TO  GRAPHYTIC  CARBON 

PRODUCTION  FROM  FIRST  CHEMICAL  AND  PHYSICAL  PRINCIPLES 
AND  FROM  EMPIRICAL  DATA  FOR  A  WIDER  RANGE  OF  FIRE  AND 
FUEL  CONDITIONS. 

E.  TO  EXAMlNINE  IN  MORE  DETAIL  THE  PHOTOCHEMICALLY  ACTIVE 

SPECIES  OF  EMISSIONS  PRODUCED  FROM  OPEN  FIRES  AND  THE 
EFFECT  OF  HIGHER  INTENSITY  COMBUSTION  CONDITIONS  ON 
THESE  EMISSIONS. 

F.  TO  MEASURE  EMISSIONS  FROM  A  WIDER  RANGE  OF  MASS  IGNITED 

PRESCRIBED  FIRES. 

G.  TO  COMPLETE  AN  EMISSIONS  CHARACTERIZATION  DATABASE  OF 

EXISTING  EMISSIONS  DATA. 


THE  FATE  OF  WILDFIRE  SMOKE 
T.  Y.  PALMER  AND  W.T.  KREISS 
Horizons  Technology,  Inc. 

San  Diego,  CA  92111-1682 

Smoke  from  the  Great  Borneo  Fire  convection  columns  behaved 
in  two  ways.  Under  stable,  dry  conditions  they  merged  and  spread 
large  distances  laterally.  When  there  was  sufficient  moisture 
for  condensation,  thev  remained  more  seperated  and  farmed 
vigorous  cumulo-form  clouds.  This  dichotomy  gives  a  basis  for 
forming  statistical-alternative  hypothesis  for  observed 
atmospheric  temperature  anomolies. 

Statistical  analvsis  shows  a  strong  heating  pulse  in  the 
troposphere  ir.  the  smokey  regions,  while  strong  cooling  occured 
in  the  lower  stratosphere.  The  alternative  postulate  for  the 
cooling  pattern  is  based  upon  past  e/periments.  Thev  indicate 
that  there  are  two  kinds  or  smoke  particles.  The  most  common  is 
a  hydrophylic  mixture  or  water,  sugars  and  other  soluables  from 
pyrolysis.  Rarer  are  the  hydrophobic^  particles,  which  are  made  up 
of  oily,  surface  act i . -  -rials  whi ch  are  probat i >  assoc i sted 
with  carbon  particles. 

The  hydrophylic  par  as  participate  in  the  usual  cloud 

processes.  Hydrophobic  .cles  "bounce"  oft  oreci pi tatx on ,  do 

not  participate  in  nuclei .  urocesses  and  probably  are  carried 

to  high  altitudes  by  cc:  .  ■ i ve  clouds.  Hvdrophylic  particles 
carried  to  these  altitude  1  lose  their  water  by  evaportaion. 

In  both  cases  a  layer  o •  tiveiv  active  sub-micron  particles 

will  be  left  in  the  strat  re. 
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Aailriii  £i  Aerosols  from  filrgfhlM  fill! 


D.  E.  Fields.  E.  T.  Arakawa,  N.  C.  Buncick,  J.  B.  Davidson,  R.  V.  Doans.  J.  Gondonnet, 
L.  D.  Hnlett.  G.  D.  Kerr,  C.  V.  Miller.  A.  M.  Soloaon,  L.  J.  Standley,* 

G.  L.  Vaughan.**  sad  M.  G.  Tslciatss 

Oak  Ridge  National  Laboratory*** 

Oak  Ridge.  TN  37831 


ABSTRACT 

Ve  hare  obtained  saaples  of  aerosols  deposited  daring  the  Hiroshima  fire 
storm  that  was  initiated  by  the  atonic  bomb  detonated  on  Angnst  6,  1943. 
These  particles,  which  we  extracted  from  streaks  of  "black  rain"  found 
on  a  plaster  wall,  are  being  studied  to  resolve  uncertainties  in  several 
areas.  Questions  being  considered  include  the  following:  What  is  the 
distribution  of  particle  sizes,  and  whet  are  the  optical  properties  of 
these  particles?  How  can  these  data  farther  our  understanding  of 
nuclear  winter  effects?  Can  examination  of  the  shape,  polymorphism,  and 
chemical  composition  of  these  particles  provide  information  about 
burning  and  particle  generation  during  fire  storms?  Can  examination  of 
residual  radionuclides  in  the  particles  and  measurement  of 

thermoluminescence  in  the  plaster  wall  yield  data  that  permit 

calculation  of  the  skin  dose  to  individuals  who  were  directly  exposed  to 
the  black  rain? 

Initial  studies  confirm  that  the  artifact  appears  to  contain  aerosol 
particles  that  we  believe  are  representative  of  the  aerosols  that  may 
lead  to  a  nuclear  winter.  Gamma  spectroscopy  measurements  indicate  the 
presence  of  natural ly-occuring  radionuclides  K-40,  Ra-226,  Ra-228,  Th- 
232,  and  Th-234,  along  with  the  fission  product,  Cs-137.  Sooty 
particles  of  varying  sizes  down  to  the  sub-micron  size  range  have  been 
detected  using  optical  photomicrography.  We  have  examined  the  particles 
using  1-ray  induced  X-ray  fluorescence  and  have  detected  Ca,  Sr.  Ba,  Fe. 
and  Zn,  with  elemental  composition  ratios  representative  of  Hiroshima 
soil.  Particle  composition  and  size  studies  using  electron  and  optical 
microscopy,  LASER- flow  cytometry.  and  electron- induced  X-ray 

fluorescence  are  continuing,  as  is  examination  of  the  artifact  using 
autoradiography  and  thermoluminescence. 
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broadcast  prescribed  bums  of  logging  slash  In  western  Washington  and  western 
Oregon  using  the  CMB  method  Include  emission  factors  for: 

1. )  Particulate  matter  by  size  class* 

a. )  PM  (4.2  to  34.2  g/kg), 

b. )  PM2.5  (1.2  to  18.2  g/kg) 

c. )  Light  absorption  properties 

2. )  Carbon  content  of  PM2.5 

a. )  elemental  carbon  (2.6S  to  40.3%) 

b. )  organic  carbon  (16. 8%  to  62. 6*) 

c. )  total  carbon  (38.7%  to  65. 2%) 

3. )  Gases  (CH-C1,  COS*  NO.*  CO*  CH.»  CO.*  and  others) 

a. )  CH.C1  (.004  to  .052  g/kg)1 2 3 4 

b. )  CCT (28.3  to  459.6  g/kg) 

c. )  CH-  (0.7  to  10.6  g/kg) 

d. )  CO4  (1014.1  to  1781.4  g/kg) 

4. )  PolynucTear  organic  materials  (BCa3P) 

a. )  BCalP  (18.8  to  398.3  5g/kg) 

b. )  Total  PAH  fraction  (0.88  to  14.11  mg/kg) 

5. )  Aldehydes*  organic  acids*  and  others. 

a. )  Acetaldehyde  (.08  to  14.13  g/kg) 

b. )  Furfuraldehyde  (.03  to  2.09  g/kg) 

We  have  proposed  making  measurements  of  most  of  the  above  emissions  for  a 
series  of  fires  (4  to  6  separate  fires)  of  chaparral  fuels  similar  to  the  fuel 
that  will  be  burned  during  the  large-scale  experiment  at  the  San  Dimas 
Experimental  Forest.  The  objectives  for  this  research  are: 

1. )  To  measure  emission  factors  for  the  above  emissions. 

2. )  To  compare  these  emission  factors  by  major  fuel  type  with  the 

results  for  the  larger  experiment. 

3. )  To  project  for  the  larger  fire  the  effect  of  the  fuel  and  fire 

variability  on  the  source  characteristics  of  emissions  sampled 
by  the  airborne  systems. 

4. )  To  test  for  similarity  of  source  characteristics  between 

airborne  and  surface  measured  emissions. 
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AIRBORNE  LIDAR  CHARACTERIZATION  OF  OPTICAL 
AND  PHYSICAL  PROPERTIES  OF  NUCLEAR  INDUCED  ATMOSPHERIC  EVENTS 


Edward  E.  Uthe 
Warren  B.  Johnson 
SRI  International 
Menlo  Park,  Ca.  9*1025 

ABSTRACT 

Experimental  data  are  needed  for  the  development  and  validation  of  climatic 
models  to  predict  the  consequence  of  atmospheric  events  triggered  by  a  nuclear 
conflict.  A  series  of  field  tests  are  planned  to  acquire  this  necessary  data 
in  the  presence  of  large  scale  smoke  and  dust  events.  We  believe  that  air¬ 
borne  lidar  observations  can  provide  information  on  the  temporal  and  spatial 
distribution  of  particulate  backscatter  and  extinction  coefficients  that  can 
not  be  obtained  by  other  means.  Path-integrated  extinction  (transmission) 
data  for  dense  aerosol  distributions  can  be  derived  from  ground-reflected 
lidar  signatures.  Smoke  and  dust  cloud  transport  and  diffusion  parameters  can 
be  derived  from  observed  geometric  distributions  of  enhanced  scattering 
volumes.  On-board  radiometric  sensors  can  supplement  the  lidar  observations 
in  evaluating  smoke  and  dust  properties,  and  the  density  of  intervening  atmos¬ 
pheric  elements  such  as  high-altitude  cirrus  clouds. 

This  paper  presents  data  examples  from  previous  studies  that  demonstrate  the 
capabilities  of  lidar,  and  proposes  a  methodology  for  the  use  of  airborne 
lidar  on  future  field  studies. 
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AIRBORNE  LIDAR  CHARACTERIZATION  OF 
OPTICAL  AND  PHYSICAL  PROPERTIES  OF 
NUCLEAR-INDUCED  ATMOSPHERIC  EVENTS 


Edward  E.  Uthe 
Warren  B.  Johnson 

SRI  International 
Menlo  Park,  CA  94025 
*415/859-4667 


This  paper  presents  data  examples  from  previous  field  experiments 
that  demonstrate  expected  capabilities  of  airborne  lidar  for 
evaluating  smoke  and  dust  source  terms  over  large  regional  areas. 
A  bibliography  of  recent  published  results  is  included. 
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SRI  INTERNATIONAL  LIDAR  SYSTEMS 


GROUND  MOBILE 

•  MARK  IX 

•  4-A. 

•  UV  DIAL 

•  IR  DIAL 

•  NEAR  IR  DIAL 

AIRBORNE 

•  ALPHA-1 

•  ALPHA-2 

•  ALARM 


0.69  |im 

0.53,  1.06,  3.8,  10.6  jim 
UV  (2) 

9.5-1 1  jim  (4) 

3-4  jim  (2) 


0.53,  1.06  jam 
0.53,  1.06  jim 
9.5-11  urn  (2) 


1.06Mm  WAVELENGTH 


20B 


CROSS-PLUME  STRUCTURE  DERIVED  FROM  ALPHA-1 
PLUME  BACKSCATTER  AT  1.06-jum 


AIRBORNE  LIDAR  OBSERVATION  OF  AN 
AEROSOL  PLUME  GENERATED  FROM 
AGRICULTURAL  BURNING  (Clear-Air  Near- 


DISTANCE  —  approximately  4  km/div 


AIRBORNE  LIDAR  MAPPING  OF  SAHARA  DUST  LAYER  OVER  CARIBBEAN  REGION 

(1969  BOMEX  PROJECT) 


•  • 


•  • 


DEVIATIONS  FROM  BEST  FIT 
EXPONENTIAL  ATMOSPHERE,  Relative  DB  Units 


I  Volcanic  dust  layer 

—  Tropopause 
"’""Cirrus  cloud  layer 


SRI  MARK  IX  LIDAR  OBSERVATIONS 
OF  MT.  ST.  HELEN'S  AEROSOL  LAYER 
Menlo  Park,  California,  2  June  1980 
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LIDAR-OBSERVED  CIRRUS  CLOUD  STRUCTURE,  MENLO  PARK,  CALIFORNIA 
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TRANSMISSION  —  percent  MAXIMUM  BACKSCATTER 


TIME - MINUTES 


CONTOUR  MAP  OF  1.06  pm  WAVELENGTH  EXTINCTION  COEFFICIENT  (km'1) 
DISTRIBUTION  ALONG  THE  LIDAR  OPTICAL  PATH  AS  A  FUNCTION  OF  SMOKE 
EVENT  TIME.  CORRESPONDING  TRANSMISSION  HISTORY  IS  PLOTTED  ABOVE 
THE  CONTOUR  MAP. 


TRANSMISSION  HISTORY  AND  OPTICAL  DEPTH  RELATIONSHIPS  FOR  SMOKE 

(ZUNI  ROCKETS) 


XTINCTION  TO  MASS  RATIO,  nl M 


TRANSMISSION  —  %  ALTITUOE  —  600  m/div  ALTITUDE 


RELATIONSHIP  OF  EXTINCTION  RATIO  AT  WAVELENGTHS 
OF  1.045-  AND  0.514-pm  AND  MEAN  PARTICLE  DIAMETER 
DERIVED  FROM  LABORATORY  RESULTS4 


F  *  FLY  ASH 

H  •  hci/nh4oh 
R  -  ROSIN  SMOKE 
P  -  SMOKE  POT 
S  •  GROUND  SILICA 
I  •  IRON  OXIDE 


/ 


o 

S  0.6 
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«  0.4 
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PARTICLE  DIAMETER.  D3J  - 

RELATIONSHIP  OF  EXTINCTION  RATIO  AT  WAVELENGTHS 
OF  1.045-  AND  0.514-*<m  AND  MEAN  PARTICLE  DIAMETER 
DERIVED  FROM  LABORATORY  RESULTS4 


21 8 


K'.v'.V V‘  O  '.' 


TRANSMIT 

1.06  pm  (IR) 

0.53  pm  (Green) 

RECEIVE 

1.06  pm  (IR) 

^  0.60  pm  (Orange) 

RESOLUTION 
3  m  vertical 
10  m  horizontal 
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An  airborne,  autotracking,  external -head  sunphotometer  for 
plume  studies. 

P. B. Russell,  T. Matsumoto,  R. F. Pueschel 
NASA  Ames  Research  Center,  Moffett  Field,  CA  94035 

An  airborne,  automatically  sun-tracking  sunphotometer 
|  consisting  of  a  detector  module,  temperatrure  control 

system,  nitrogen  purge  system,  mechanical  drive  chain  and 
data  collection  system  has  been  developed.  The  instrument 
measures  continuously  vertical  profiles  of  optical  depth  and 
transmissivity!  The  detector  head  is  mounted  outside  the 
aircraft  cabin  to  maximize  sun-viewing  opportunities  and  to 
avoid  viewing  through  a  cabin  window  at  varying  transmission 
angles.  Using  6  separate  detectors,  each  with  its  own 
optical  filter,  rather  than  a  single  detector  with  a  filter 
wheel,  measurements  at  6  different  wavelengths  from  0. 38  to 
1. 02  urn  can  be  taken  simultaneously  to  investigate  plumes  of 
even  limited  extent  at  aircraft  speeds  as  high  as  100  m/sec. 

The  data  collection  system  consists  of  a  multiplexer,  12-bit 
analog- to-digital  converter,  aircraft  data  interface 
electronics,  and  a  Hewlett-Packard  HP9816  computer  with  a 
floppy  disk  and  a  printer.  The  sampling  period  is 
synchronized  to  the  aircraft  data  system  that  updates  the 
navigation  data  approximately  every  2  seconds.  The  6 
detector  signals,  detector  temperature,  pressure  altitude, 
longitude,  latitude,  tracking  error,  sun  tracker  azimuth  and 
elevation,  and  GMT  are  recorded  on  floppy  disks  and  also 
printed  out  for  data  backup.  The  HP9816  is  also  used  for 
postflight  data  processing  and  Langley  and  optical  depth 
plots. 

Mounting  configurations  for  Convair  C990  and  C131, 

MacDonnel -Douglas  DCS  and  Lockheed  C130  aircrafts  have  been 
developed  or  are  being  contemplated. 
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Biologically  Important  Issues  for  Research  on  the 
Physical  Effects  of  Nuclear  War 

Marie  A.  Harwell 
Cornell  University 


*  Endpoints  of  Concern  -  Human  Effects 

*  Fallacy  in  Delaying  Biological  Analyses 

*  Uncertainties  -  Biological  vs.  Physical  Issues 

*  Vulnerabilities  Define  Relevance  of  Research 
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Key  Conclusions  from  SCOPE-ENUWAR  Study 


Indirect  Effects  >  Direct  Effects 
Ecological  Effects 

-  Unprecedented  in  Extent  and  Intensity 

-  Different  Ecosystems  with  Different  Vulnerabilities 

-  Long-Term  Recovery  Prospects 

Natural  Ecosystems  -  Human  Carrying  Capacity  <  1  %  Global  Population 
Agroecosystems  Most  Vulnerable  to  Disturbance 
Agricultural  Effects 

-  Temperature  Most  Important 

-  Few  degrees  vice  tens  of  degrees 

-  Many  mechanisms  for  temperature  effects 

-  Other  Climatic  Effects  Also  Important 

-  If  No  Climatic  Effects,  Still  Problem  Human  Subsidies 
Possibility  of  Elimination  or  Severe  Disruption  in  Grain  Production 

-  for  at  least  N.Hemisphere 

-  for  at  least  one  growing  season 
Food  Limitations  Critical 

-  Stores  Depleted  Prior  Next  Growing  Season 

-  Stores  Could  Support  Only  Small  Fraction  of  Global  Population  for  1  Yr 
Other  Effects 

-  Radiation 

-  global  fallout  -  little  effect 

-  local  fallout  -  large  effect 

-  UV-B  -  potentially  large  effect 

-  Air  Pollutants  -  potential  local  effects 

-  Synergisms 

Societal  Effects  -  little  studied,  potentially  large  effects 


Mechanisms  for  Climatic  Effects  on  Agriculture 


*  episodes  of  freezing  during  growing  season 

-  acute  and  chronic  period 

-  variance  vs.  average  temperature 

-  spatial  extent 

-  long-term  reverberations 

-  poorly  estimated 

*  shortening  of  growing  season 

-  chronic  period  -  historical  records 

-  variance  vs.  average  temperatures  -  statistical  analyses 

-  spatial  extent  -  latitudinal  gradients 

-  poorly  estimated  re  nuclear  war 

*  increase  in  growing  season  length  requirements  by  crop 

*  insufficient  thermal  time 

-  chronic  period 

-  average  only 

-  better  estimates 

*  insufficient  integrated  insolation  time 

-  chronic  period 

-  average  only 

-  daylength 

-  little  attention 

*  suboptimal  environmental  conditions 

-  temperature 
-light 

-  precipitation  **  key  issue 


Variance  Issues 


*  Acute  Period 

-  spatial  heterogeneity 

-  temporal  heterogeneity 

-  extremes  vs.  averages 

*  Chronic  Period 

-  extremes  vs.  averages 

-  variable  vs.  fixed  growing  season 

-  diurnal  range 


Issues  re  Averages  in  Temperatures  Over  Time 


*  Uniform 

1. )  reduce  daily  averages 

-  unknown  effect  re  diurnal  range 

2. )  reduce  daily  max,  mins 

-  no  change  in  diurnal  range 

3. )  reduce  daily  min,  not  max 

-  increase  diurnal  range 

4. )  reduce  daily  max,  not  min 

-  decrease  diurnal  range 

♦Non-Uniform 

1 . )  increase  frequency  of  cold  events 

2. )  increase  duration  of  cold  events 

3. )  increase  intensity  of  cold  events 
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Research  Needs  Climatic  Effects 

Acute  Period 

*  Rate  of  onset  of  temperature  and  light  reductions 

*  Absolute  temperatures  vice  decrements 

*  Spatial  extent  of  freezing  events 

*  Spatial  and  temporal  heterogeneity 

*  Smaller  scales  -  global/regional/local 

*  Daylength 

*  Light  spectra 

*  Precipitation 

*  Seasonality  of  effects 

*  Variance  issues 


Research  Needs  Climatic  Effects 


Chronic  Period 


*  Duration,  intensity  of  temperature  decreases 

-  frost-free  season  duration 

-  effects  of  seasonality 

-  Southern  Hemisphere  effects 

-  surface  water  effects 

*  Duration,  intensity  of  light  decreases 

-  daylength 

-  spectral  quality 

*  Duration,  intentity  of  precipitation  decreases 

-  coastal,  erogenic  effects 

-  monsoon  disruption 

*  Variance  issues 

-  spatial  heterogeneity 

-  average  changes  vs.  variance  changes 

-  uniform  vs.  non-uniform 

-  diurnal  ranges 

*  Ocean/atmospheric  interactions 

-  coastal  vs.  continental  effects 

-  ocean  currents 

*  Long-term  feedbacks 

-  albedo 

-  Antarctic  ice  pack  dynamics 

-  greenhouse  gases 

-  biological  feedbacks  (e.g.,  fires,  CO2,  hydrology, 

desertification) 
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Research  Needs  Other  Effects 

Radiation 

*  Distribution  and  levels  of  local  fallout 

-  variety  of  scenarios 

-  spatial  heterogeneity 

-  increased  emphasis  Beta-emitters 

*  Internal  doses 

-  fate  and  transport 

-  routes  to  humans  post-war  conditions 

-  dose  levels 

-  gamma.  Beta,  and  alpha-emitters 

UV-B 

*  Increased  levels  expected 

-  duration 

-  spatial  distribution  -  latitude/altitude 

-  effects  particulates  on  ozone  layer 


Air  Pollutants 


NOx 

acid  precipitation 

hydrocarbons 

others 


*  Variety  of  pollutants 

-CO 

-  pyrotoxins 

-  ozone 

-  asbestos 

*  Levels  expected 

-  spatial  distribution,  extent 

-  duration 


L  I^A'1 


I  » 


SAIC /MCLEAN,  VA 


SOVIET  RESEARCH  ON  GLOBAL  EFFECTS 
(PAPER  GIVEN  AT  THE  DNA  GLOBAL  EFFECTS  PROGRAM  TECHNICAL  MEETING, 

AMES  RESEARCH  CENTER,  25-27  FEBRUARY  1986) 

There  are  obvious  reasons  for  interest  in  Soviet  research  into  the  global 
effects  and  consequences  of  the  nuclear  war. 

First,  for  some  time  now  Soviet  scientists  have  been  active  participants 
in  the  international  scientific  research  efforts  to  study  the  global  effects 
and  consequences  of  a  nuclear  war.  Indeed,  it  could  be  said  that  this  has 
been  a  major  area  of  at  least  apparent  Soviet  scientific  cooperation  with 
the  international  scientific  community.  Soviet  scientists  have  also  appeared 
as  actors  in  the  public  debate  in  the  Free  World  concerning  the  consequences 
and  potential  implications  of  the  effects.  While  the  participation  of  Soviet 
scientists  in  these  studies  has  been  generally  welcomed  by  the  international 
scientific  community  for  scientific  as  well  as  political  reasons,  it  is 
useful  to  assess  the  actual  contribution  which  Soviet  scientists  have  made 
in  this  field  and  the  validity  of  their  claims  to  have  independently  validated 
the  key  findings  and  projection  of  their  western  colleagues. 

Second,  as  we  all  know,  many  issues  concerning  global  effects  are  fraught 
with  large  uncertainties  and  are  quite  sensitive  to  assumptions,  choices  of 
scenarios,  value  judgements,  etc.  The  Soviet  Union  has  many  highly  qualified, 
indeed  outstanding  scientists  who,  one  might  expect,  could  make  significant 
contributions  from  an  independent,  non-Western  perspective  on  the  basis  of 
their  own  data  and  experiments  to  our  understanding  of  the  effects  and  their 
consequences . 

Third,  in  principle,  findings  on  the  global  effects  could  have  important 
implications  for  the  nuclear  powers'  defense  policies,  force  level  requirements 
views  on  deterrence  and  strategic  stability,  targeting  doctrines,  weapons 
development  programs,  arms  control,  and  so  on,  or  in  brief  on  the  feasibility 
and  risks  of  resort  to  nuclear  weapons.  It  is  important,  therefore,  to  try  to 
ascertain  actual  Soviet  views  on  and  degrees  of  real  Soviet  concern  about  the 
various  possible  global  effects  and  their  implications  for  Soviet  war  fighting 
doctrine  and  strategy.  In  this  respect  a  distinction  may  have  to  be  made 
between  what  the  Soviets  regard  as  serious  problems  and  those  which  they  seek 
to  exploit  for  policy  purposes .  It  is  also  obvious  that  asymmetries  in  East- 
West  understanding  and  evaluations  of  the  consequences  of  nuclear  war  could  be 
dangerously  destabilizing. 

Fourth,  it  is  inevitable  that  perceptions  of  or  claims  about  various 
global  effects  and  their  consequences  will  be  used  for  propaganda  and  political 
purposes.  It  is  especially  important  to  keep  in  mind  that  in  the  Soviet  Union 
science  and  scientists  are  under  the  control  of  the  state  and  that  scientists 
are  called  upon  by  the  party  leadership  to  use  their  influence  and  prestige  in 
support  of  the  Soviet  Union's  propaganda  and  foreign  policy  objectives.  This 
was  clearly  spelled  out  by  Boris  Ponomarev  at  the  all-union  conference  of 
Soviet  scientists  held  in  Moscow  in  May  1983.  One  result  of  this  is  that  the 
better  known  or  publicized  Soviet  work  on  global  effects  in  the  West  has  been 
in  a  large  measure  deliberately  targeted  on  Western  audiences  not  just  for 
scientific  but  also  for  political  effects.  Indeed,  a  significant  part  of 
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this  work  is  not  published  in  the  Soviet  Union  or  appears  only  in  publications 
specifically  aimed  at  foreign  audiences.  This  raises  questions  about  the 
seriousness  of  Soviet  scientific  purposes  in  this  field  and  Soviet  objectives. 
At  the  same  time,  however,  there  is  other  Soviet  work  which  is  given  less 
publicity  abroad,  but  which  may  reflect  more  serious  Soviet  interests  and 
concerns  and  which,  therefore,  may  deserve  more  careful  attention  on  our  part. 

Let  me  say,  that  as  a  "Sovietologist"  or  political  scientist  and  not  a 
physical  scientist,  I  am  only  an  observer  of  and  not  a  direct  participant  in 
the  scientific  effort  to  study  global  effects  of  nuclear  war.  I  have,  however, 
a  long-standing  interest  in  studying  Soviet  views  on  nuclear  war  fighting 
and  its  consequences,  which  I  have  sought  to  derive  from  studies  of  Soviet 
military  doctrine  and  strategy,  targeting  and  weapon  employment  concepts 
and  especially  of  the  Soviet  civil  defense  program,  which  deals  with  Soviet 
views  on  the  effects  of  a  nuclear  war  and  the  character  of  the  post-attack 
environment . 

My  interest  in  the  Soviet  approach  and  contribution  to  the  so  called 
Nuclear  Winter  issue  began  with  my  participation  at  the  1983  International 
Seminar  on  Nuclear  War  held  in  Erice,  Italy,  at  which  V.  Aleksandrov  of  the 
Computing  Center  of  the  USSR  Academy  of  Sciences  first  presented  the  findings 
of  the  first  Soviet  model  of  the  climatic  consequences  of  nuclear  war.  It 
was  further  spurred  by  the  conversations  I  had  there  with  Aleksandrov  and  the 
Vice-President  of  the  USSR  Academy  of  Sciences,  Ye.  Velikhov.  This  led  in 
1984  to  a  quick  survey  sponsored  by  DNA  on  Soviet  studies  and  exploitation  of 
the  Nuclear  Winter  hypothesis  covering  the  years  1983-1984.  The  primary 
objective  of  my  study  was  to  examine  the  question  whether,  as  the  Soviets 
publicly  claimed,  Soviet  scientists  had  in  fact  independently  verified  the 
findings  and  predictions  of  the  Crutzen-Birks  paper  published  in  1982  in 
Amblo  and  more  importantly  those  of  the  so  called  TTAPS  study,  published  in 
1983. 


•Vn 

as 


& 


>>■ 

$ 

1 


rv  1 


i 


My  conclusions  were  that  there  was  no  real  basis  for  this  Soviet  claim.  What 
the  Aleksandrov  and  Stenchikov  1983  computer  model  had  done  was  to  uncritically 
use  the  scenarios  and  assumptions  about  amounts  of  fire,  soot  and  dust  and 
their  distribution  from  the  Amblo  and  TTAPS  studies.  Aleksandrov  acknowledged 
that  he  had  obtained  the  TTAPS  findings  in  April  1983  during  a  presentation 
made  of  them  in  Cambridge,  Massachuettes ,  and  had,  simply  used  the  TTAPS 
assumptions  not  for  the  base-line  case  but  for  a  10,000  MT  exchange  scenario. 

The  initial  findings  were  further  distorted  by  an  error  of  a  factor  of  more 
than  three  in  interpreting  the  radiative  effects  of  the  TTAPS  estimate  of 
optical  depth  by  treating  dust  as  smoke  and  by  assuming  that  smoke  would 
absorb  100  percent  of  the  incident  solar  radiation,  which  was  the  equivalent 
of  assuming  an  injection  of  some  1,000  Tg  of  smoke  Instead  of  the  225  Tg 
assumed  in  the  TTAPS  study. 
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A  great  deal  has  been  written  about  the  shortcomings  and  errors  of  the 
original  Aleksandrov/Stenchikov  computer  simulation  model  and  the  assumptions 
underlying  it.  The  Soviets  were  slow  in  acknowledging  and  correcting  the 
errors  in  their  computations  and  findings.  Subsequently  Aleksandrov  did 
modify  his  model,  at  least  to  some  extent.  He  did  look  at  the  100  MT  TTAPS 
case,  at  moving  smoke  and  scavenging,  the  transport  of  smoke  to  the  southern 
hemisphere,  and  reduced  somewhat  the  amount  of  injected  smoke  in  his  com¬ 
putations.  The  model  also  assumed  significant  heating  at  higher  elevations, 
such  as  the  Himalayas  and  other  mountain  peaks.  Overall,  it  could  be  said 
that  the  Soviet  contribution  to  the  1983-1984  Nuclear  Winter  research  remained 
limited  to  the  Aleksandrov/Stenchikov  computer  simulation  model  using  TTAPS 
scenarios  and  assumptions,  which  represented  a  seriously  flawed  effort,  It 
was  quickly  outdistanced  by  modeling  done  in  the  West .  The  Aleksandrov/ 
Stenchikov  work  was  initially  only  published  in  English,  and  to  this  day 
continues  to  be  cited  by  its  English  title  in  Soviet  Russian  language  pub¬ 
lications.  An  abridged  version  of  this  model  was  eventually  published  in 
1984  in  a  Soviet  mathematics  journal.  As  to  the  Aleksandrov/Stenchikov  1984 
paper  "Update  on  the  Climatic  Impacts  of  Nuclear  Exchange,"  it  seems  to  have 
been  generally  ignored  in  Soviet  publications  and  articles  on  Nuclear  Winter, 
and  appears  to  have  had  no  significant  effect  on  Aleksandrov's  presentation 
at  a  meeting  of  East-European  scientists  which  was  held  in  November  1984. 

There  is  little  doubt  that  during  1983-1984  the  Soviets  sought  to  gain 
maximum  propaganda  value  from  their  more  extreme  findings  on  Nuclear  Winter 
effects,  just  as  they  did  from  other  studies  of  global  medical  and  biological 
effects  on  nuclear  war  using  worst-case  Western  scenarios  (WHO).  The  basic 
message  which  Soviet  propaganda  sought  to  convey  was  that  nuclear  war,  and 
especially  Nuclear  Winter  is  not  survivable  and  that  so-called  limited  use 
of  nuclear  weapons  (a  concept  attributed  to  the  U.S.)  would  lead  to  as  serious 
catastrophic  consequences  as  a  full  scale  nuclear  exchange.  Without  ques¬ 
tioning  the  sincerity  of  individual  Soviet  scientists,  one  could  come  to  the 
conclusion  that  from  the  Soviet  viewpoint  the  Nuclear  Winter  issue  was  treated 
less  as  a  serious  scientific  issue  requiring  careful  scientific  inquiry, 
although  one  must  recognize  the  limitations  of  the  computational  resources 
at  the  Academy's  Computing  Center,  than  as  a  propaganda  opportunity  in  support 
of  the  Soviet  peace  campaign. 

It  may  be  worth  noting  that  in  1984  some  Soviet  academicians  writing  in 
domestic  Soviet  publications,  expressed  "indignation"  at  what  they  claimed  to 
be  attempts  by  some  Western  scientists  to  moderate  the  worst-case  predictions 
of  the  Nuclear  Winter  phenonemon .  They  condemned  these  attempts  "to  cast 
doubts  on  the  conclusions  reached"  by  the  "joint  efforts"  of  Soviet  and 
American  scientists  who  were  said  to  have  "firmly  established  the  existance 
of  such  global  threshold  phenomena."  The  articles  protested  the  alledged 
efforts  of  some  Western  scientist  to  "discredit  the  results  of  the  work  of 
conscientious  researchers  (presumably  including  Aleksandrov/Stenchikov)  who 
are  opening  people’s  eyes  to  the  danger  threatening  them."  /For  example,  see 
Academician  V.  Goldansky  and  Professor  S.  Kapitsa,  "To  Prevent  a  Catastrophe," 
Izvest lya,  25  July  1984^7 
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Although  on  various  occasions  a  number  of  Soviet  scientist  members  of 
various  Sovie.  institutes  have  expressed  interest  in  the  Nuclear  Winter  issue, 
in  practice  there  was  little  forthcoming  from  them  during  1983-1984.  Indeed, 
despite  some  Soviet  improvements  in  their  modeling  they  did  not  eliminate  some 
of  their  serious  errors  until  1985.  After  a  year  or  so  of  so-called  Soviet 
collaboration  with  the  International  Nuclear  Winter  Research  effort.  Dr.  Turko 
noted  that,  aside  from  Aleksandrov's  computational  modeling  work,  there  had 
been  "no  substantial  physical  data  and  little  evidence  of  objective  scientific 
analysis  forthcoming"  from  the  Soviet  side.  A  public  request  by  Dr.  Sagan 
in  November  1983  to  Soviet  scientists  for  data  on  the  particle  size  distri¬ 
bution  function  of  debris  from  Soviet  atmospheric  nuclear  weapon  tests  before 
1963,  for  information  on  particle  size  and  absorbtion  coefficients  from  large 
fires  in  the  Soviet  Union,  and  for  Soviet  views  on  the  more  likely  nuclear 
war  scenarios  has  remained  essentially  without  response.  True  in  1984  and 
1985  there  has  been  talk  of  the  Soviets  conducting  various  fire  experiments, 
but  there  is  no  indication  that  these  have  been  carried  out.  Indeed  in  the 
summer  of  1984  the  Soviets  indicated  that  the  modeling  of  large  fires  was 
beyond  the  capacity  of  the  computers  at  the  Academy's  Computing  Center.  The 
Soviets  say  that  they  have  some  data  on  fires  but  have  difficulties  inter¬ 
preting  them.  Nevertheless  one  might  note  that  the  Soviet  Nuclear  Winter 
studies  make  no  reference  to  or  use  of  the  data  on  the  massive  Siberian 
forest  fires  in  the  summer  of  1915  which  are  discussed  in  Soviet  fire  litera¬ 
ture.  No  mention  is  made  of  other  large  forest  and  peat  fires,  such  as 
those  in  1972.  In  both  cases,  very  large  areas,  (2.6  million  square  miles 
in  1915)  were  subjected  to  heavy  persistent  smoke. 

Soviet  research  into  and  publications  on  global  effects  offer  a  mix  of 
serious  scientific  inquiries  and  what  may  be  more  properly  characterized  as 
science  propaganda.  The  public  aspects  of  the  Soviet  work  do  not  reflect  the 
existence  of  private  differences  of  opinion  among  Soviet  scientists  and  skep¬ 
ticism  about  some  of  the  published  findings,  which,  however,  become  evident 
in  non-public  or  private  discussions  between  American  and  Soviet  scientists. 

Of  course,  no  public  debate  of  the  sort  practiced  in  the  West  is  permitted  in 
the  Soviet  Union,  especially  on  issues  which  Soviet  propaganda  likes  to  exploit. 

Soviet  work  on  global  effects,  and  especially  on  Nuclear  Winter,  suffers, 
at  least  as  far  as  its  public  aspects  are  concerned,  from  a  number  of  con¬ 
straints.  One  of  these  is  the  absence  of  Soviet  war  scenarios  and  the  pro¬ 
hibition  on  public  discussion  of  potential  effects  of  nuclear  strikes  on  the 
Soviet  Union  in  terms  of  potential  targets  which  may  be  attacked,  likely  urban 
and  forest  fires,  fallout  patterns  and  so  on.  True,  Soviet  literature  on 
weapon  effects  and  civil  defense  discusses  fires  and  fire  propagation,  fallout , 
radioactive  contamination,  average  fuel  loading  in  cities  and  fuel  loading 
requirements  for  fire  storms  and  so  on,  but  these  are  discussed  in  the  abstract 
and  never  in  terms  of  specific  Soviet  cities, targets  or  regions.  When  it  comes 
to  such  specifics,  the  Soviets  fall  back  on  citing  Western  data  on  what 
could  happen  to  the  U.S.  or  Western  Europe  in  the  event  of  an  attack.  Even 
when  the  Soviets  illustrate  the  point  that  a  2  C  decline  in  ambiant  temperature 
will  damage  crops,  they  cite  the  case  of  such  an  effect  on  the  Canadian  wheat 
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crop.  Furthermore,  for  propandanda  reasons  the  Soviets  tend  to  rely  on  worst 
case  Western  scenarios,  no  matter  how  unrealistic.  It  is  difficult  to  believe, 
however,  that  there  are  no  Soviet  scenarios  and  independent  assessments  of  the 
effects  of  nuclear  strikes  in  the  light  of  such  scenarios.  Unfortunately, 
these  are  not  shown  to  the  West,  and  indeed,  they  are  likely  to  be  classified. 

Another  constraint,  or  so  we  are  led  to  believe,  is  the  limitation  on 
the  capabilities  of  Soviet  computers,  which,  however,  does  not  explain  the 
tenacity  with  which  the  Soviets  have  perpetuated  some  of  the  errors  in  their 
computations  over  a  number  of  years.  Still  another  factor  appears  to  be 
bureaucratic,  i.e.,  the  question  of  who  is  going  to  lead  the  research  effort, 
which  institute  will  be  given  a  leading  role,  and  how  this  research  fits  into 
the  planning  of  the  work  of  the  Soviet  science  establishment.  During  1983- 
1984  the  leading  role  was  played  by  a  small  group  of  people  in  the  Academy's 
Computing  Center,  under  the  leadership  of  V.  Aleksandrov,  who,  were  only 
interested  in  the  mathematical  modeling  of  atmospheric  and  biological  pro¬ 
cesses,  but  had  little  expertise  in  atmospheric  physics,  fire  phenomenology, 
chemistry,  etc.  and,  therefore,  were  only  too  willing  to  use  uncritically  Weste 
data.  The  involvement  of  other  institutes  with  staffs  having  such  qualifica¬ 
tions,  such  as  the  Academy's  Institute  for  Atmospheric  Physics  or  Dr.  Yuri 
Israel's  Laboratory  for  Maintaining  of  Natural  Environment  and  climate, 
to  data  has  been  limited.  While  Velikhov  may  have  been  directing  the  overall 
effort,  the  bureaucratic  problem  appears  to  have  become  seriously  aggravated 
by  V.  Aleksandrov's  mysterious  disappearance  in  Madrid  in  March  1985  during  a 
speaking  engagement,  given  that  he  was  the  leader  of  the  main  Soviet  Nuclear 
Winter  effort  and  its  primary  spokesman  in  the  international  scientific 
community.  It  is  noteworthy,  that  in  the  most  recent  Soviet  English  language 
publication  on  Nuclear  Winter,  The  Night  After,  Aleksandrov's  name  has  dis¬ 
appeared  from  all  references  to  his  1983  study.  In  effect,  poor  Aleksandrov 
may  be  becoming  a  non-person  in  his  own  homeland,  which  may  also  suggest  that 
international  collaboration  in  Nuclear  Winter  research  may,  under  some  condi¬ 
tions,  not  be  devoid  of  serious  hazards  to  the  health  of  Soviet  scientists. 

It  has  been  pointed  out  to  me,  however,  that  in  a  recently  published  monograph 
by  Stenchikov  and  Dr.  Carl  of  the  GDR  Academy  of  Sciences,  issued  by  the  latter 
Aleksandrov's  name  was  cited  in  the  references.  We  may  be  seeing  a  case  of  his 
resurrection,  at  least  in  name,  in  the  GDR  if  not  the  USSR. 

It  should  be  emphasized  that  Soviet  scientists  are  not  newcomers  to 
research  into  the  effects  of  nuclear  war.  They  have  carried  out  such  work 
for  a  long  time.  Their  work,  however,  has  tended  to  focus  on  the  distribution 
of  radioactive  fallout,  the  problems  of  long  term  contamination  by  radio¬ 
isotopes,  the  biological  and  ecological  effects  of  radiation,  and  the  effects 
of  the  introduction  into  the  atmosphere  of  various  chemical  aerosols,  etc. 

A  reading  of  Soviet  publications  also  reflects  a  strong  concern,  still  ongoing, 
in  the  problem  of  depletion  of  the  ozone  layer.  Soviet  scientists  have  brought 
up  the  point  about  sunlight  absorbtion  by  nitrogen  oxides,  said  to  have  been 
observed  during  Soviet  atmospheric  nuclear  tests  during  the  1960s,  and  about 
the  effect  of  large  amounts  of  greenhouse  gases  released  into  the  lower  atmos¬ 
phere,  which,  it  is  said,  could  lead  to  localized  "Nuclear  Summer"  effects. 

It  may  well  be  that  the  Soviets  have  more  to  contribute  in  the  field  of 
atmospheric  chemistry  and  effects  of  injections  of  various  aerosols  then  in 
the  modeling  of  Nuclear  Winter. 
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While  we  continue  to  see  a  mix  of  Soviet  science  research  and  propanganda, 
there  is  some  recent  evidence  of  Soviet  attempts  to  provide  more  realistic 
stimulations,  and  to  rely  on  less  extreme  scenarios,  although  these  continue 
to  be  borrowed  from  Western  sources.  For  example,  in  the  book.  The  Night  After, 
published  in  1985,  we  find  a  slick  propaganda  piece  by  Velikhov  along  with  a 
revised  model  by  Stenchikov  which  he  did  jointly  with  Carl.  Among  other 
changes  the  Stenchikov  computations  use  an  initial  absorbtion  optical  depth 
of  2.2  rather  than  6  as  in  the  old  model,  more  scavenging,  narrower  initial 
latitudinal  spread  of  aerosols  and  so  on,  which  also  results  in  a  reduction 
in  the  calculated  temperature  changes,  and  brings  them  more  in  accord  with 
the  Western  studies.  In  a  chapter  by  Yu  Israel,  he  uses  a  base  line  scenario 
of  5,000  MT  and  notes  a  low  probability  of  the  injection  of  significant  amounts 
of  aerosol  particles  from  forest  fires  into  the  stratosphere.  At  the  same 
time,  however,  in  a  press  conference  given  by  Stenchikov’s  chief.  Academician 
N.  Moiseyev  in  November  1985  he  kept  citing  the  findings  of  the  old  Soviet 
models  rather  than  those  of  Stenchikov.  Also  in  1985  a  pamphlet  on  Disastrous 
Effects  of  Nuclear  War  -  Soclo-Economical  Aspects  was  published  primarily  in 
English,  French,  German  and  Spanish,  under  the  auspices  of  the  Soviet  Scientific 
Research  Council  on  Peace  and  Disarmament  with  an  editorial  board  headed  by 
Academician  P.N.  Fedoseyev.  The  authors'  claim  about  the  socio-economic 
consequences  of  nuclear  war  are  clearly  predetermined  by  their  claim  that 
"the  main  target  in  a  nuclear  war  are  people  and  their  civilization,  not  the 
military  systems  of  the  states". 

The  Soviet  contribution  to  the  study  of  tha  biological  effects  of  nuclear 
war  has  also  been  characterized  by  a  mix  of  science  and  propaganda.  The  latter 
has  been  particularly  noticeable  in  the  case  of  the  work  carried  out  at  the 
Computing  Center  under  the  leadership  of  Yuri  Svirezhev  and  his  attempt  to 
have  his  paper  adopted  by  SCOPE  despite  its  obvious  lack  of  realism  and  the  use 
of  obsolete  Nuclear  Winter  projections  and  other  seriously  wrong  assumptions  and 
computations  all  of  which  were  intended  to  make  the  worst  possible  case  for 
the  global  effects  of  nuclear  war. 

It  is  not  clear  at  present  how  the  Soviet  effort  to  study  global  effects 
is  organized  and  whether  it  has  a  sense  of  direction.  Although  there  are 
indications  of  interest  in  these  problems  among  Soviet  scientists,  the  Soviet 
effort  appears  to  be  rather  modest  in  scale,  and  by  no  means  all  the  work 
done  has  been  readily  accessible  to  Western  scholars.  In  the  absence  of  an 
indepth,  comprehensive  analysis,  which  remains  to  be  undertaken,  one  gains 
the  impression  that  in  terms  of  its  scope,  level  of  effort  and  richness  the 
Soviet  iNuclear  Winter  research  program  lags  far  behind  comparable  programs 
and  efforts  carried  out  in  the  West.  The  Soviet  contribution  is  still  largely 
limited  to  computer  simulations  despite  the  weaknesses  of  Soviet  computer 
capabilities.  Naturally,  for  scientific  and  all  the  more  political  reasons, 
Soviet  scientists  can  be  expected  to  continue  to  participate  in  the  inter¬ 
national  cooperative  effort  to  study  the  global  effects  and  consequences  of 
nuclear  war.  Of  course,  there  can  be  departures  from  this  practice,  such 
as  the  last  minute  failure  of  Soviet  scientists  to  attend  the  International 
Seminar  on  Nuclear  War  in  Erice  in  1985. 
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For  all  the  Soviet  public  promotion,  largely  in  the  West,  of  what  could 
be  characterized  as  a  worst  case  Image  of  the  consequences  of  a  nuclear  war, 
there  are  no  clear  answers  as  yet  to  the  questions  I  mentioned  at  the  beginning 
concerning  the  reasons  why  it  is  Important  to  assess  Soviet  research  into 
global  effects  and  to  seek  to  identify  actual  Soviet  perceptions  and  concerns 
regarding  their  implications.  My  present  intentions  are  not  only  to  update 
and  improve  my  study  for  DNA  on  Soviet  research  and  exploitation  of  Nuclear 
Winter,  but  also  to  carry  out  an  in-depth  inquiry  into  Soviet  research  and 
actual  views  on  all  aspects  of  the  post-nuclear  attack  environment .  For 
this  purpose  I  intend  to  analyze,  with  help  from  physical  scientists,  not 
only  Soviet  research  results  more  widely  circulated  in  the  West,  but  also 
Soviet  research  which  is  only  reported  in  Soviet  Russian  language  sources 
and  which  may  not  necessarily  always  be  specifically  identified  as  dealing 
with  global  nuclear  war  effects  although  it  may  be  relevant  to  them.  As  I 
have  mentioned,  there  is  a  substantial  body  of  Soviet  fire  and  nuclear  effects 
literature  and  relevant  scientific  papers  published  by  the  Academy,  the 
universities  or  appearing  in  a  variety  of  Soviet  Technical  Journals.  I  also 
intend  to  use  in  my  analysis  whatever  I  can  learn  about  what  has  transpired 
at  meetings  and  conversations  between  Soviet  and  Western  scientists.  The 
latter  information  may  be  especially  valuable  for  setting  Soviet  research 
and  views  of  the  Soviet  science  establishment  in  better  perspective  and  for 
gaining  a  clearer  appreciation  of  true  Soviet  concerns  and  priorities. 

In  my  work  so  far  I  already  owe  a  great  deal  to  the  kind  help  of  a  number 
of  scientists  who  have  not  only  provided  me  with  copies  of  Soviet  research 
papers  they  have  received  from  Soviet  scientists  but  also  notes  on  their 
conversations  with  Soviet  scientists.  I  hope  that  with  the  indulgence  of 
our  scientific  community  my  research  will  also  continue  in  the  future  to 
benefit  from  this  invaluable  assistance  to  my  further  work  in  this  field. 
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Radiological  Doae  Assessments 
From  Nuclear  War 


Charles  Shapiro 

San  Francisco  Stats  University  and 
Lawrence  Livermore  National  Laboratory 


Ted  Harvey 
Kendall  Peterson 

Lawrence  Livermore  National  Laboratory 


Abstract  of  presentation  at  the  Defense  Nuclear  Agency  (DNA)  Global  Effects  Program 
Technical  Meeting:  February  25-27,  1080;  NASA  Ames  Research  Center. 

Recently  ws  undertook  an  assessment  of  the  radiological  doss  following  a  nuclear  war 
as  part  of  the  ICSU-SCOPE-ENU WAR  study  published  as  “Environmental  Consequences 
of  Nuclear  War;  Vol.  I,  Physical  and  Atmospheric  Effects.*  The  following  are  suggested 
as  follow-on  research  topics  to  this  subject. 

1.  Evaluate  internal  dose  contributions  in  a  post-nuclear  war  environment.  For  a  normal 
atmosphere,  various  researchers  have  provided  means  to  calculate  organ  dose  for  many 
nuclides  and  food  pathways.  These  calculations  will  not  apply  to  an  environment 
perturbed  by  “nuclear  winter*  effects. 

2.  Extend  the  calculations  of  global  fallout  to  include  the  stratospheric  contribution  in  a 
smoke-perturbed  atmosphere  using  a  fully  interactive  GCM. 

3.  Improvements  to  the  local  fallout  code  KDFOC2  and  its  documentation. 

4.  Improve  understanding  of  the  radiological  doae  commitments  associated  with  the  po¬ 
tential  targeting  and  damaging  of  civilian  and  military  nuclear  fuel  cycle  facilities. 

5.  GLODEP2  code  improvements  and  scenario  calculations  for  global  fallout.  Code  im¬ 
provements,  parameter  sensitivities  and  scenario  studies  will  be  conducted. 

6.  “Reconcile”  local  fallout  code  differences.  To  investigate  and  communicate  on  local 
fallout  code  differences  and  establish  benchmark  type  problems  related  to  nuclear  win¬ 
ter. 

This  work  was  performed  under  the  auspices  of  the  U.S.  Department  of  Energy  by  the 
Lawrence  Livermore  National  Laboratory  under  Contract  W-7405-Eng-48. 
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MODELS 


GL0DEP2:  Estimates  gamma  ray  dose 

from  global  fallout. 
Authors:  L.  Edwards,  T.  Harvey, 

K.  Peterson 

Reference:  UCID-20033  (1984) 


XFU^rpnif  iyyi 

S 


GRANTOUR:  Calculates  the  3-D  transport, 

diffusion,  and  wet  and  dry  £ 

deposition  and  dose  from  tropo-  £ 

spheric  global  fallout. 

Authors:  J.  Walton  and  M.  MacCracken 
Reference:  UCID- 19985  (1984)  gjj 

KDF0C2:  Calculates  local  fallout. 

Authors:  T.  Harvey  and  F.  Serduke 
Reference:  UCRL-52858  (1979) 


MAIN  POINTS 


LOCAL  FALLOUT  CAN  COVER  LARGE  AREAS  OF  PROTAGONIST  COUNTRIES 
WITH  LETHAL  LEVELS  OF  RADIATION  (LARGE  UNCERTAINTIES). 

GLOBAL  FALLOUT  IS  MORE  HAZARDOUS  THAN  PREVIOUSLY  THOUGHT 
(•-  3-*. 3  0  RADS  IN  N.H.  MID  LATITUDES) 

N.H.  AVERAGE  IS  ABOUT  15  RADS 
S.H.  AVERAGE  IS  ABOUT  0.5  RADS 

A  "NUCLEAR  WINTER"  ATMOSPhERE  REDUCES  GLOBAL  FALLOUT  DOSES 
BY  15°:  IN  N.H.,  AND  INCREASES  THE  DOSE  MARGINALLY  IN  S.H. 

(to  <  1  RAD  )  . 

"HOTSPOTS"  MO6  KM2)  ARE  ABOUT  6  'R  7  TIMES  HIGHER  THAN 
THE  N.H.  AVERAGES. 

SMALLER  YIELD  DEVICES  RELATIVELY  LARGER  DOSES. 

T  E  GLOBAL  FALLOUT  POPULATION  DOSE  IS  ABOUT  € X 1  0 1  0 
PERSON-RADS. 

A  "WORST  CASE"  SCENARIO  ATTACK  ON  NUCLEAR  FUEL  CYC.  r 
FACILITIES  WOULD  GREATLY  EXACERBATE  THE  FALLOUT  F- 

WE  HAVE  TREATED  MAINLY  ExiERNAL  GAMMA  DOSE. 

INTERNAL  DOSES  FROM  GAMMAS,  EETAS  AND  ALPHAS  APE  A 
SIGNIFICANT  BUT  MORE  DIFFICULT  TO  ESTIMATE.  BETA- 
CAN  HAVE  SIGNIFICANT  IMPACT  ON  BIOTA. 
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Scenario  A  Scenario  B 

\nox  (1983)  5300  Mt  TTAPS  (Turco  et  al.  1983)  5 

baseline  nuclear  war  reference  nuclear  war 


surface  explosions 
surface  explosions 


TABLE  3.7.  Global  fallout  dose  assessments  (rads)  (unperturbed  atmosphere  no 


8  l 

a  o 

• — '  3 

*  s 

u  * 
<0  S 

*  u 

U  <0 
=3  « 

c  73 

<u  3 
C  c 

1  8 

^  5 

-Q  £; 
-  ,2 
:>  2 


„  NCaNO’ftdHt'.^H  .r 

oa  ®  r:  ^  • — *  o  o  o  o  p  d  ^  .r 


conneno^iflNXNoo 
w*  ^HNN-HobodNcid 

^  N  N  — 1 


n  rtssiOrtN'fel-iNMa 

cq  NNiiddoddddoN 

CM  M  — t 


^  as’reo’feiflNHifloo 

< *  N-isdoddddddd 

^  IM  CM  — i 


_  t»5osn®i|'ftiq^'foo 

m  cioo-iddddddoida 
n  •»  — • 

_  U)rtOI®00(OflOlflHNOf 

•fi^Nddodddddoo 

N  M  —i 


n  pi  -r 

a.  0,  O 
tj  id  H 
Q  Q 
O  O  < 

»— 3  Ch 

U  U  o 

M  M  50 


.  Jj  « 

-r* 

7  I  2  o 
■  /•  -  ~o  f> 
7.  /)  v  -o 


bO 

bO 

c 

c 

C 

C 

*7 

’  7 

•7 

7 

3 

3 

3 

3 

C 

C 

c 

c 

.2 

3 

O 

o 

c 

w 

w 

•-» 

«-* 

*-» 

U 

U 

V 

u 

V* 

•  n 

• 

•  — > 
G 

•  n 

*  ** • 

■  "* 

•  -* 

■  — ■ 

• 

u 

u 

u 

t- 

u 

a> 

V 

G 

z 

2  <5  Z  o  -/5  •/:  t.  d 

Q000-*0000 
yi  t'  i/5  n  m  n  ci 
I  I  I  I  V)  i  l  i  i 
000030033 
r^toco  — •  — *  — *  root— 


G  3  3 

rr>  n  r>  G  .r  CJ 


<  a  <  O  < 


WV.AV.A' 


Figure  7.5  Global  fallout:  ac-  ted  whole  body  gamma  dose  (rad.)  from 
6235  explosions  totaling  2031  Mi.  fission  product*  (scenario  A).  An  8  day 
tropospheric  deposition  decay  constant,  characteristic  of  a  winter  injection, 
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TABLE  8.9.  Global  fallout  dose  using  three-dimensional  CRANTOUR  (summer 
scenario).  Comparison  of  nuclear  winter  atmosphere  (smoke)  and  unperturbed  at¬ 
mosphere  (no  smoke).  Doses  are  in  rads.  Because  GRANTOUR  only  calculates 


uox,  1983),  unpert 
nox,  1983),  perturl 


Dose  (rads) 


Figure  7.6  Comparison  of  radionuclide  global  dose  distribution  for  eases  with  unper¬ 
turbed  and  smoke-perturbed  climates  (tropospheric  contributions  only). 


LOCAL  FALLOUT 


Table  7.1  Surface-burst  warheads  in  a  phased  nuclear  exchange.  All  weapons 
are  assumed  to  be  50%  fission  yield. 


Weapon 
yield  (Mt) 

Number  of  warheads 

Initial 

counter¬ 

force 

phase 

Extended 

counter¬ 

force 

phase 

Industrial 

counter- 

value 

phase 

Final 

phase 

Full 

baseline 

exchange 

0.05 

0 

300 

0 

250 

550 

0.1 

975 

150 

50 

8 

1183 

0.2 

0 

250 

50 

121 

421 

0.3 

500 

250 

0 

125 

875 

0.5 

1000 

200 

0 

25 

1225 

1.0 

250 

495 

160 

125 

1030 

5.0 

0 

50 

15 

8 

73 

Total  surface- 

burst  yield 

-1000 

-1000 

-250 

-250 

-2500 

*»»  M 


Figure  7.1  48-hour  dose  predictions  for  a  1-Mt  all-fission  weapon  detonated 
at  the  surface.  A  mid-continental  Northern  Hemisphere  summer  wind  profile 
was  used.  The  double-lobed  pattern  is  due  to  a  strong  directional  wind 
shear  that  is  typical  during  this  season.  For  a  1-Mt  weapon,  the  lofting  of 
radioactivity  is  so  high  that  topographic  features  are  not  expected  to  play  a 
large  role  in  pattern  development;  thus,  a  fiat  surface  has  been  used.  The 
protection  factor  is  1.  The  local  terrain  is  assumed  to  be  a  rolling  grassy 
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Initial 

counter- 

force 

phase 

Extended 

counter¬ 

force 

phase 

Industrial 

counter- 

value 

phase 

Final 

phase 

Full 

baseline 

exchange 

Europe 

Eastern 

0 

2.9 

0.6 

0.8 

4.3 

U.S.S.R. 

Western 

0.5 

0.5 

0.1 

0.2 

1.3 

U.S.S.R. 

1.6 

2.3 

0.7 

1.7 

6.3 

Eastern  U.S. 

0 

4.7 

1.0 

1.4 

7.1 

Western  U.S. 

4.4 

2.3 

0.7 

6.6 

8.0 

fi*"'*  9  8  Gamma-ray  dose  rate  area  integral  versus  time  after  shutdown  or  detonatk 
(Chester  and  Chester.  1976). 


1-MT  weapon  Reprocessing 
plant  10-year 

- — high  level  storage 


30-Day  high 
level  storage- 


3600-MW  (th)  LWR 

(steady  state)-^T 

10- Year  storage 
at  reactor 

Gamma-ray  dose  rate  1  meter  above  a 
smooth  plane  for 

(1)  1-MT  fission  weapon 

(2)  1200-MWfe)  LWR  power  reactor 

(3)  10- Year  storage  of  spent  fuel  at  the 
reactor 

1500-Metric  ton/year  fuel  reprocessing 
plant  with 

(4)  10-Year  high  level  storage 

(5)  30-Day  high  level  storage 

vs 

Time  after  shutdown  or  detonation 


>-  3 


~  0J 

c  a 
a  o 


Time  (sec 


quarter  years 


•  Forage-cow-milk, 

•  Pasture-  meat  (including  steer,  lamb,  and  hog 
meat), 

•  Leafy  vegetable, 

•  Grains  and  other  non-leafy  crops, 

•  Fresh  and  salt  water  fish,  and  other  edible 
aquatic  animals. 


Surface  air  concentration  of  the  radionuclide, 
Deposition  of  Cs-137  on  the  grass  or  dry  feed, 
Fraction  of  the  nuclide  that  remains  on  the 
forage. 

Residence  time  of  the  nuclide  on  the  forage, 
Intake  rates  of  both  wet  (grass)  and  dry  (hay) 
feed  consumed  by  the  cow, 

Concentration  in  milk. 

Dose  per  unit  activity  ingested, 

Infant  s  milk  consumption  rate. 


•  (a)  little  or  no  effect, 

•  (b)  slight  to  moderate  effect, 

•  (c)  moderate  to  large  effect, 

•  (d)  large  effect. 
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NOT  FOR  QUOTATION  OR  CITATION 

THE  SIBERIAN  MEGAFIRE  OF  1915:  A  NATURAL 
SMOKE  INJECTION  IN  THE  SAME  RANGE 
AS  "NUCLEAR  WINTER"  MODELS 

The  reduction  in  sunlight  and  temperature  over  large 
areas  presumed  to  result  from  the  20  to  650  teragrams^  of 
smoke  that  could  be  produced  by  a  major  nuclear  exchange  has 
been  extensively  modeled. 

Natural  analogs  of  real  use  in  calibrating  the  models 
have  not  been  easy  to  find.2'^ 

However#  a  uniquely  relevant  event  in  the  recent  past 
has  emerged  from  obscurity.  The  vegetation  of  central 
Siberia  is  rich  and  complex.  In  the  western  portion# 
coniferous  forest  is  interspersed  with  steppe#  birch  groves# 
and  peat  bogs  containing  a  long-term  accumulation  of 
biomass.  In  the  east#  larch  is  the  dominant  forest  species# 
and  permafrost  is  more  extensive.  In  1915#  a  prolonged 
drought#  the  worst  in  two  generations#  turned  the  Siberian 
landscape  into  tinder.  In  July  and  August#  often  violent 
wildfires  raged  from  the  Ob  to  the  Lena  (Fig.  1).  In  all# 
seven  hundred  thousand  square  miles  (roughly  1.8  million 
square  kilometers)  burned.  The  complete  devastation  of 
tracts  of  forest  as  large  as  Germany  were  reported — the 
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approximately  132*000  square  mile  belt  between  the  Angara 
and  lower  Tunguska  Rivers*  and  55*000  square  miles  of 
accessible  economic  timber*  the  most  valuable  of  which  was 
coniferous*  were  lost.  A  pall  of  smoke  overspread  an  area 
of  2.6  million  square  miles  extending  widely  beyond  the 
limits  of  the  fires*  reaching  8°  west  and  21°  east  of  the 
69°  to  112°  east  longitude  belt  in  which  burning  occurred — 
in  a  generally  northeasterly  drift.  In  the  densest  smoke* 
day  turned  to  night  and  the  cows  came  home.  Over  an  area  of 
702*000  square  miles*  visibility  dropped  to  between  14  and 
70  feet. 

Data  collected  from  350  properly  completed 
questionnaires  out  of  500  dispatched  throughout  Siberia  by 
A.  Vosnesensky  of  the  Irkutsk  magnetic-meteorological 
observatory  were  reduced  in  1916-17  by  him*  J.  Belyaeff*  and 
V.  B.  Shostakovitch*  and  a  summary  was  finally  published  in 
English*  after  a  delay  due  to  war  and  revolution*  in  April 
1925  in  the  Journal  of  Forestry.4 

While  these  smoke  emissions  were  observed  in  the  lower 
atmosphere  their  optical  depth*  scale*  spread*  and 
duration — an  average  of  51  days  —  render  their  behavior  and 
effect  on  climate*  regional  and  global*  as  interesting  as 
the  hypothetical  high-altitude  smoke  injections  employed  by 
the  several  generations  of  existing  "Nuclear  Winter"  models. 

Using  Small  and  Bush's  methodology*^  the  lower  bound  of 
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the  range  yields  about  40  teragrams  of  smoke  from  the  forest 
and  taiga  fires  reported  by  the  Shostakovitch — if  the 
reported  combustion  of  peat#  which  covers  hundreds  of 
thousands  of  square  kilometers  of  the  region#  and  humus# 
which  is  reported  to  have  burned  down  to  its  sandy  substrate 
in  some  areas#  is  completely  excluded  from  the  calculations. 

However#  using  the  methodology  employed  by  CrQtzen  and 
Birks6  or  Turco#  Toon#  Ackerman#  Pollack#  and  Sagan#7  the 
same  calculation  yields  approximately  180  to  720  teragrams 
of  smoke.  Including  peat  and  humus  could  add  an  additional 
1#000  teragrams#  for  burning  up  to  a  depth  of  two  meters# 
corresponding  to  a  fuel  loading  of  100  g/cm2  or  more  was 
reported  in  many  places.  The  50-teragram  "threshold"  for 
global  effects  may  have  been  exceeded  in  less  than  a  week  at 
the  height  of  the  fires  if  the  high  wildfire  biomass 
loadings  (.5  to  2.2  G/cm2)  they  cite  are  correct. 

"Nuclear  Winter"  effects  are  calculated  to  be  most 
severe  in  continental  interiors  in  mid-summer#^  the  season 
when  solar-heating-induced  buoyancy^  is  most  likely  to 
transport  smoke  generated  in  the  boundary  layer  into  the 
high  atmosphere  and  beyond  the  effects  of  low-altitude 
removal  processes.  In  1915#  the  dilution  of  optical  depth 
dependent  climatic  effects  was  limited  by  the  prolonged 
confinement  of  the  smoke  within  Siberia--an  area  more  than 
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an  order  of  magnitude  smaller  than  the  30°  to  60°N.  Lat. 
circumglobal  belt  considered  in  the  model  of  Turco  et  al. 
Each  100  TG  of  smoke  over  Siberia  corresponds  to  more  than 
1*000  TG  in  their  1-D  model.  This  concentration  compensates 
for  the  more  gradual  rate  of  injection.  Given  these  "worst- 
case"  aspects  of  the  real  situation  in  Siberia*  one  might 
expect  to  have  seen  climatic  and  ecological  changes  on  a 
regional  to  global  scale:^  effects  that  could  be 
significant  if  emissions  from  these  often  intense  fires  were 
on  the  level  predicted  by  Small  and  Bush's  methodology*  and 
catastrophic  if  the  model  parameterizations  of  Turco  et  al.* 
the  NAS*  or  SCOPE  are  in  good  correspondence  with  the 
vicissitudes  of  natural  history.  For*  regardless  of  the 
altitude  at  which  the  modelers  elect  to  inject  the  smoke 
they  model*  it  too  must  have  its  source  on  the  ground  and 
its  origin  in  the  boundary  layer. 

It  is  significant  that  data  from  Siberian 
meteorological  stations  record  the  coexistence  of  widespread 
smoke  and  sporadic  darkness  with  a  hot*  dry  summer; 
insolation  was  reduced  by  an  average  of  15%  in  July  and  35% 
in  August.  These  observations  are  consistent  with  the  early 
and  unpublished  results  of  3-D  modeling  by  MacCracken  et 
al.*  in  which  standard  aerosol  injection  parameters — 
altitudes  below  1.5  km  above  surface  were  applied  to 
"Nuclear  Winter"  aerosols  (the  model  predicted  surface 


S' 159 (5)  -  REVISED  DRAFT  SUBMITTED  TO  NATURE  (2/23/86) 


warming — personal  communication/  M.  MacCracken).  While  this 
was  sufficient  to  retard  the  ripening  of  cereal  crops  by  10 
to  15  days/  despite  the  expectation  of  an  early  harvest  on 
account  of  the  elevated  temperatures  and  low  rainfall/  it 
scarcely  coincides  with  the  collapse  of  the  agricultural 
ecosystem  postulated  by  Harwell's  model*!  of  the  response  of 
biota  to  comparable  transients  in  the  optical  depth  of  the 
atmosphere. 

Much  needs  to  be  done  to  better  quantify  the 
conflagration  of  1915.  Should  Vosnesensky's  data  base/  which 
was  stored  at  Irkutsk  circa  1925/  prove  inaccessible/  a 
wealth  of  dendrochronological  evidence  still  survives  in  the 
trunks  of  mature  trees  in  Siberia  and  could  be  used  to 
augment  the  written  record  of  the  events  in  question.  The 
eyewitnesses  who  remain  are  nonagenarians  whose  oral  history 
must  be  gathered  soon  or  not  at  all. 

Until  such  data  are  forthcoming/  one  of  the  most 
edifying  experiments  available  to  the  climate  modeling 
community  might  consist  in  loading  the  3-D  interactive 
models*^  with  boundary  layer  smoke  of  optical  depth 
corresponding  to  the  1915  reported  observations  over  the 
appropriate  areas  of  Siberia.  One  could  then  allow  the 
computers  to  run  on  historical  emissions  and  meteorology*3 
rather  than  purely  parametric  injections.  For  the 
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mathematically  perfect  homogeneity  of  the  smoke  in  the  3-D 
global  climate  model  grids  is  very  different  from  the 
fractal  complexity  of  clouds  in  the  real  world. ^  will  the 
models  depict  the  collapse  of  the  Hadley  Circulation15  and  a 
frigid  end  for  the  tropical  rainforests  and  ricefields  of 
the  Third  World?  Or  will  1915r  as  modeled*  resemble  the 
year  itself?  A  slightly  warmer  than  average  year  in  which 
the  echo  of  the  Guns  of  August  drowned  out  the  cries  of  fire 
in  uncrowded  Siberia  and  the  world's  attention  turned  to  the 
terrible  realities  of  the  Great  War. 

The  author  wishes  to  thank  E.  Bauer*  G.  Carrier* 

M.  Harwell*  A.  Hecht*  and  R.  P.  Turco  for  their  stimulating 
discussions*  M.  Kuchement  of  the  Harvard  Russian  Research 
Center  for  his  aid  in  clarifying  the  geographic  extent  of 
the  fires*  and  C.  Sagan  for  first  proposing  the  utility  of 
fire  data  from  the  Soviet  Union  in  "Nuclear  Winter" 
research. 15 
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SURFACE  TEMPERATURE  EFFECTS  OF  FOREST  FIRE  SMOKE  PLUMES 


A  proposal 


by 

Alan  Robock 

Department  of  Meteorology 
University  of  Maryland 
College  Park,  Maryland  20742 


Although  it  is  impossible  to  study  an  actual  nuclear 
winter,  nature  provides  us  with  analogs,  occurrences  in  the 
actual  climate  system  that  can  teach  us  about  some  of  the 
postulated  interactions.  Almost  every  year  recently,  large 
forest  fires  have  burned  in  North  America,  Asia  or  the 
tropics  that  have  produced  extensive  smoke  plumes.  Some  of 
these  have  been  measured  to  have  optical  depths  in  the 
visible  larger  than  3,  as  thick  as  postulated  nuclear  winter 
smoke  clouds.  Although  nuclear  winter  smoke  plumes  would  be 
a  mixture  of  forest  fire  smoke  and  smoke  from  burning  cities 
(in  some  scenarios)  and  hence  have  different  optical 
properties,  we  can  learn  much  that  is  relevant  to  nuclear 
winter  by  the  surface  air  temperature  responses  to  forest 
fire  smoke  plumes.  The  Great  Smoke  Pall  of  1950,  resulting 
from  Canadian  forest  fires,  was  reported  to  have  lowered 
daytime  surface  air  temperatures  in  Washington,  D.C.,  by  6 
degrees  C  after  atmospheric  residence  of  one  week.  Detailed 
studies  of  contemporary  plumes  with  adequate  data  sets  will 
allow  us  to  attempt  to  verify  such  effects. 

It  is  proposed  here  to  analyze  surface  air  temperature 
response  to  the  presence  of  smoke  plumes  in  the  atmosphere. 
Cases  will  include  forest  fire  plumes  with  sources  in 
Western  Canada,  Montana,  Siberia,  New  Guinea  and  Borneo  that 
have  been  identified  on  satellite  imagery  during  the  past  5 
years.  Surface  temperatures  under  the  plume  will  be 
compared  to  those  outside  me  plume  by  various  techniques, 
including  the  succesful  me  mod  of  using  MOS  errors  that 
detected  surface  temperatum  fleets  of  the  Mount  St.  Helens 
dust  cloud.  Results  will  r  ••  compared  to  calculations  from 
radiation  models. 
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Abstract 

A  key  limitation  of  standard  numerical  models  of  nuclear  winter  is  the 
indirect  pararnetrisation  of  the  all  important  (soot)  cloud-radiation 
interaction,  and  the  associated  assumption  of  homogeneous  clouds  at  sub- 
grid  scales  However,  basic  ideas  in  turbulence  suggest  that  soot  wiU  be 
concentrated  by  cascade  type  processes  into  highly  inhomogeneous  patchy 
regions  characterised  by  strong  fluctuations  of  the  form  Pri&C.ac)  •.£c‘a 
Tu!  the  probability  distribution  of  a  soot  concentration  f luctuation(AC) 
exceeding  a  fixed  threshold  _l  Using  various  tracers,  we  empirically 
3fe  ’f  e  exponent  v  ■  3  whereas  for  temperatures  in  normal 
climate  xv.T  -  5  Padiation  balance  arguments  suggest  that  in  nuclear 

winter  fluctuations  in  soot  are  directly  related  to  surface  temperature 
fluctuations  and  tnat  nuclear  winter  may  lower  a  a,t  from  ^  to  3  yie'dirg 

much  more  intermittent  temperature-1  ■  vve  argue  4  hat  tor  Mo  logic  a1 
systems  the  intermittenry  of  the  temperature  may  be  as  important  the 
mean  smre  extreme  low  temperatures-  even  if  only  bnefly  at  tamed-  can  dn 
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great  damage.  We  conclude  that  the  patchiness  of  nuclear  winters  could 
make  them  more  severe  than  previously  believed. 


1. Introduction: 

Due  to  their  limited  range  of  scales.  General  Circulation  Models  are  not  able 
to  handle  directly  the  central  physical  processes  in  Nuclear  Winter  (Crutzen 
and  Birks.  1982,  Turco  et  al.  1983,  Alexandrov  arid  Stenchikov.  1934,  Covey 
et  al.,  1984).  In  particular,  the  all-important,  soot  clouds  are  highly 
parametrised,  and  consequently  the  soot-radiation-ternperature 
interactions  are  modelled  very  indirectly.  This  is  very  unfortunate,  since 
the  biological  consequences  of  a  Nuclear  Winter  are  expected  to  be  highly 
dependent,  on  the  intensity,  duration,  and  frequency  of  occurence  of  extreme 
low  temperature  episodes-  not  only  on  the  average  temperature  drop 
(Harwell  and  Hutchison  1986)  The  fundamental  significance  of  such  extreme 
variability  have  been  widely  discussed  by  biologists,  under  the  rubric 
"variance  effects".  Below,  we  directly  address  this  question  by  relating 
temperature  variability  to  the  patchiness  of  the  soot  clouds 

2  How  patchiness  is  produced 

The  patchiness  of  clouds  can  be  understood  m  the  context  of  turbulent 
cascades  In  (scaling)  cascade  processes  a  given  quantity  'the  flux  of 
energy  for  the  dynamics,  here  the  soot)  is  concentrated  step  ty  M.ep  Corn 
to  the  smallest  scale  by  a  similar  process  af  each  scale  The  era' i  ef ?e •  • 
of  such  processes  is  to  concentrate  quantities  in  /ery  s  af  w  ■  ■!  -  :  *  *  *  - 


Fig  !  give  such  an  example  a  cen  -of  sire  >'n  -  ’s  split,  due  to  non-linear 

interactions  with  cells  of  roughly  the  same  size  in  the  sub-cells  -  of  size 
/n4l  =  >. l n ,  l  oei rig  the  '.constant .■  scale  ratio  of  the  discretised  cascade- 

arid  its  content  is  randomly  distributed  among  the  different,  sub-cells,  in 
such  a  way  as  to  conserve  (on  average)  the  cascade  quantities  at  each  scale. 
In  the  homogeneous  case  (Fig!  a),  the  quantity  is  evenly  distributed  among 
the  different  sub-cells  and  no  variations  in  concentration  occurs  However, 
with  a  simple  modification  (Fig  i  b)  "dead/alive"  sub-cells  (no  soot/  soot) 
strong  patchiness  is  produced,  arid  the  number  of  the  "alive”  sub-cells  of 
size  l  varies  as: 

N(  e )  H>s,  Dj<  d 

d  being  the  dimension  of  the  space  on  which  the  cascade  occurs  (e  g  3,  or  4 
in  space-time  space),  and  Ds  is  the  fractal/Hausdorff  dimension  of  the 

support  of  "activity",  j.e.  where  soot  is  present  and  even  so  dense  that  its 
density  is  singular  in  respect  to  the  usual  volume  (Lebesgue)  measure 

As  soon  as  we  abandon  the  alternative  "dead/alive"  (corresponding  t. :  ♦-■e  , 

model,  Frisch  et  a 1 . ,  1978)  to  the  more  realistic  alternative  wv  ■  • 
(the  "a-model",  Schertzer  and  Lovejoy,  1983.  1968a  b;  a  .»>:’»  *  •  -  * 

(decreasing)  dimensions  of  the  (sparser  arid  ;;  a*'  ••• 
singular  behaviour  of  the  density  of  the  soot 

3  E  x  t  r  e  m  e  flu  c  t  u  a  1 1  o  r:  a :  /  .  ^ ■ 
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statistical  moments  of  the  soot  occurs.  The  physical  reason  is  that  the  the 
density  of  the  soot  is  so  singular  that  its  integration  over  low-dirnensionnal 
sets  is  insufficient  to  smooth  out  this  singular  behaviour,  of  the  high  order 
statistical  moments.  This  type  of  divergence  is  associated  with  algebraic 
fall-off  of  probability  distributions,  we  may  therefore  assume  that  the  the 
probability  distribution  of  the  soot  concentration  (or  of  its  fluctuations) 
will  have  "fat  tails",  e.g.: 

Pr(AC(Ax)>Ac)  ^  Ac-^;  Ac»1. 

where  AC(Ax)  are  fluctuations  of  the  concentration  C  of  the  soot  on  a 
distance  Ax,  and  Ac's  are  (high)  thresholds,  a  is  the  intermittency  exponent 
(independant  of  the  spatial  or  temporal  scale).  One  may  note  that  many 
meteorological  fields  have  been  investigated  in  this  manner,  the  rain-field 
(Love joy,  1981,  Lovejoy  and  Schert2er1985),  dynamics  and  potential 
temperature  (Schertzer  and  Lovejoy,  1983,  1985a)  and  surface  temperature 
(Ladoy  et  al.,  1985,  Lovejoy  and  Schertzer,  1986,  Lovejoy  et  al.,  1986). 
Furthermore,  even  in  strongly  anisotropic  frameworks,  a  wide  variety  of 
continuous  cascade  processes  are  possible,  (Schertzer  and  Lovejoy,  1985b, 
1986. 

As  appropriate  soot  dispersion  data  were  not  available,  analysis  of  aerosol 
tracers  were  used  to  characterise  tracer  cloud  patchiness.  Below,  data 
from  SF6  and  CBr2F2  tracer  experiments  (Wilson  et  al.  1976)  were  analysed 

(Visvanathan,  1985,  Schertzer  et  al.,  1986).  Measuring  stations  were  placed 
along  circles  of  radii  100m,  200m,  400m,  centered  on  a  point  source  (see 
Fig. 2).  The  scaling  and  intermittent  behaviour  of  the  concentration  of  these 
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tracers  could  therefore  be  studied  over  a  range  of  several  hundred  meters. 
Though  the  amount  of  data  was  limited,  the  intermittency  exponent  of  the 
fluctuations  of  the  concentration  was  estimated  as: 

&Ac  ^  3 

in  support  of  the  scaling  hypothesis  this  asymptotic  behaviour  of  the 

probability  distributions  was  apparently  the  same  for  the  different 

distances  studied  (fig.  3).  The  assumption  of  scaling  behaviour  over  a  much 

wider  range  of  time  and  space  scales  is  justified  by  the  observations  of 

Doury,  1980  that  the  r.m.s.  width  of  clouds  using  a  wide  variety  of  tracers 

(including  nuclear  clouds),  obey  algebraic  laws.  This  means  that  it  is  not 
unreasonable  to  expect  <x£c  to  be  constant  over  a  wide  range  of  scales. 

However,  large  scale  fire  experiments  will  be  necessary  to  fully  confirm 
this  idea. 


In  the  radiative  transfer  equation,  if  we  consider  the  emissivity 
fluctuations  to  be  of  second  order,  with  respect  to  fluctuations  of 
concentration,  we  obtain  a  simple  relation  between  the  fluctuations  of 
surface  temperatures  (AT)  and  soot  concentration  (AC): 

AT/T  ^  -AC/4. 

This  crude  approximation  is  consistent  with  different  radiative  transfer 

models  (e.g.  MacCracken  1983,  Penner  and  Haselman  1983),  and  implies  that 
the  intermittency  exponent  of  the  surface  temperature  (a^T)  is  the  same  as 

that  of  the  concentration  (a^c).  In  normal  conditions  (Ladoy  et  al ,  1985, 
Lovejoy  and  Schertzer,  1986,  Lovejoy  et  al.,  1986)  we  find  empirically: 

Abnormal  ^  ^ 


we  therefore  expect  far  more  interrmttency  in  Nuclear  Winter  conditions, 
since  the  present  study  suggests: 

^Nuclear  Winter^  3 

This  represents  a  much  more  frequent  occurence  of  extreme  temperatures. 


in  the  near  future  more  sophisticated  studies  will  undoubtly  be  performed 
by  coupling  radiative  transfer  schemes  with  cascade  models  of  clouds  (see 
for  example  Gabriel  et  al„  1986). 


We  have  briefly  indicated  how  turbulent  cascade  models  provide  a 
framework  in  which  the  extreme  fluctuations  of  a  Nuclear  Winter  could  be 
studied  both  theoretically  and  experimentally.  In  particular,  they  allow  us 
to  directly  relate  the  intermittent  behaviour  of  the  surface  temperatures  to 
the  patchiness  of  soot  clouds. 


Using  SF6  and  C8r2F2  as  tracers,  a  preliminary  data  analysis  enabled  us  to 
estimate  the  value  of  the  intermittency  exponent  of  the  fluctuations  of  the 
concentrations  of  the  tracers.  Using  a  crude  approximation  of  the  radiative 
transfer  equation,  we  argued  that  under  nuclear  winter  conditions,  where 
temperature  fluctuations  are  likely  to  be  controlled  by  fluctuations  in  the 

radiation  field  (and  this  in  turn  by  clould  patchiness),  that  this  exponent 
should  be  roughly  the  same  as  that  of  the  surface  temperatures  ((<XAT)NuclMr 

winter  'v  3)  Since  *he  comparable  value  for  normal  conditions  is  (ct^-i-)^.^,  ^ 
5,  we  have  preliminary  quantitative  evidence  that  Nuclear  Winter  will  be 
characterised  by  more  frequent  occurences  of  extreme  temperature 


fluctuations.  Large  scale  controlled  fires  will  be  indlspenslble  in  testing 
these  conclusions. 


The  biological  consequences  of  a  Nuclear  Winter  are  likely  to  depend 
prlmarly  on  the  occurences  of  extreme  conditions,  and  not  so  much  on  the 
average  drop.  We  therefore  conclude  that  Nuclear  winter  is  likely  to  be 
more  severe  than  previously  thought. 
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Figure  1 :  A  schematic  representation  of  how  turbulence  cascades 
treat  the  break-up  of  a  single  eddy  (represented  by  the  central 
square)  via  non-linear  interactions  during  a  single  step  in  the 
cascade  process.  Both  schemes  shown  here  are  isotropic,  the  left 
hand  side  is  homogenous,  and  the  right  hand  side  intermittent. 
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AIR  FORCE  OPERATIONAL  TEST  &  EVAL  CTR 
ATTN:  MAJ  N  RODRIGUES 

AIR  FORCE  SPACE  DIVISION 

ATTN:  YNC  CAPT  K  O’BRYAN 
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AIR  FORCE  TECHNICAL  APPLICATIONS  CTR 
ATTN:  J MARSHALL 

AIR  FORCE  WEAPONS  LABORATORY,  NTAAB 
ATTN:  CAPTLEONG 
ATTN:  JJANNI 
ATTN:  J  W  AUBREY,  NTED 
ATTN:  LT  LAHTI 
ATTN.  LTCOLV  BLISS 

AIR  UNIVERSITY 

ATTN:  LTCOL  F  REULE 

AIR  WEATHER  SERVICE,  MAC 

ATTN:  MAJ  J  SCHLEHER 

BALLISTIC  MISSILE  OFFICE/DAA 
ATTN:  LT  ROTHCHILD 
ATTN:  MYSP/CAP  TOMASZEWSKI 

DEPUTY  CHIEF  OF  STAFF/XOX 
ATTN:  AFXOX 

STRATEGIC  AIR  COMMAND 

ATTN:  CAP  CONNERY 

STRATEGIC  AIR  COMMAND/XPXF 
ATTN:  T  BAZZOLI 

DEPARTMENT  OF  ENERGY 

ARGONNE  NATIONAL  LABORATORY 
ATTN:  H  DRUCKER 
ATTN:  M  WESLEY 
ATTN:  P  BECKERMAN 

BROOKHAVEN  NATIONAL  LABORATORY 
ATTN:  B  MANOWITZ 
ATTN.  E  WEINSTOCK 

DEPARTMENT  OF  ENERGY 
ATTN:  INEDDOW 
ATTN:  T  HARRIS 

DESERT  RESEARCH  INSTITUTE 
ATTN:  JHALLETT 
ATTN:  J  HUDSON 

LAWRENCE  BERKELEY  NATIONAL  LAB 
ATTN:  H ROSEN 

UNIVERSITY  OF  CALIFORNIA 

LAWRENCE  LIVERMORE  NATIONAL  LAB 
ATTN:  CRMOLENKAMP 
ATTN:  C  SHAPIRO 
ATTN:  F LUTHER 
ATTN:  GBING 
ATTN:  G  SIMONSON 
ATTN:  JPENNER 
ATTN:  J  POTTER 
ATTN:  L-10  A  GROSSMAN 


ATTN:  L-262  A  BROYLES 
ATTN:  L-262  J  KNOX 
ATTN:  L-453  L  ANSPAUGH 
ATTN:  M  MACCRACKEN 
ATTN:  R  MALONE 
ATTN:  RPERRET 
ATTN:  SGHAN 

LOS  ALAMOS  NATIONAL  LABORATORY 
ATTN:  DSAPPENFIELD 
ATTN:  E  J  CHAPYAK 
ATTN:  E JONES 
ATTN:  E  SYMBALISTY 
ATTN:  GGLATZMAIER 
ATTN:  GM  SMITH 
ATTN:  L  HAUER 
ATTN:  L  CLOUTMAN 
ATTN:  PHUGES 
ATTN:  T  YAM  ATT A 

OAK  RIDGE  NATIONAL  LABORATORY 
ATTN:  D  FIELDS 

SANDIA  NATIONAL  LABORATORIES 
ATTN:  A  L  JOHNSON 
ATTN:  BZAK 
ATTN:  DDAHLGREN 
ATTN:  D  FORDHAM 
ATTN:  D  WILLIAMS 
ATTN.  KD  BERGERON 
ATTN:  LTROST 
ATTN:  MD  BENNETT 
ATTN.  R  C  BACKSTROM 

OTHER  GOVERNMENT 

CENTRAL  INTELLIGENCE  AGENCY 
ATTN:  AWARSHAWSKY 
ATTN:  R  NELSON 

DEPARTMENT  OF  AGRICULTURE 
ATTN:  D  HAINES 

DEPARTMENT  OF  TRANSPORTATION 
ATTN:  COLMROESCH 

ENVIRONMENTAL  PROTECTION  AGENCY 
ATTN:  R  COTHERN 
ATTN:  WE  FALLON 

FEDERAL  EMERGENCY  MANAGEMENT  AGENCY 
ATTN:  BW  BLANCHARD 
ATTN.  D  BENSON  NP-CP-MR 
ATTN:  DKYBAL 
ATTN:  J  POWERS 
ATTN:  J  RUMBARGER 
ATTN:  S  ALTMAN 
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GENERAL  ACCOUNTING  OFFICE 
ATTN:  A  PIERCE 
ATTN:  PJBOLLEA 
ATTN:  V  BIELECKI 


ATTN:  N  CRAYBILL 
ATTN:  WRCOFER 


ATTN:  RHABERLE 
ATTN:  OTOON 
ATTN:  R  YOUNG 
ATTN:  T  ACKERMAN 

NATIONAL  BUREAU  OF  STANDARDS 
ATTN:  GMULHOLLAND 
ATTN:  R  LEVINE 
ATTN:  RREHM 
ATTN:  R  SCHRACK 

NATIONAL  BUREAU  OF  STANDARDS 
ATTN:  H  BAUM 

NATIONAL  CENTER  ATMOSPHERIC  RESEARCH 
ATTN:  JKIEHL 
ATTN:  S  SCHNEIDER 
ATTN:  S  THOMPSON 
ATTN:  VRAMASWAMY 

NATIONAL  CLIMATE  PROGRAM  OFFICE 
ATTN:  AHECHT 
ATTN:  MYERG 

NATIONAL  OCEANIC  &  ATMOSPHERIC  ADMIN 
ATTN:  FFEHSENFELD 
ATTN:  JDELUISI 
ATTN:  R  DICKINSON 
ATTN:  R  PUESCHEL 
ATTN:  V  DERR 

NATIONAL  OCEANIC  &  ATMOSPHERIC  ADMIN 
ATTN:  B  HICKS 

NATIONAL  RESEARCH  COUNCIL 
ATTN:  KBEHR 
ATTN:  R  DEFRIES 

NATIONAL  SCIENCE  FOUNDATION 
ATTN:  B  BEASLEY 
ATTN:  EBIERLY 
ATTN:  HVIRJI 
ATTN:  LHAMATY 
ATTN:  R  SINCLAIR 
ATTN:  R  TAYLOR 
ATTN:  SKEENY 

NUCLEAR  REGULATORY  COMMISSION 
ATTN:  R  ALEXANDER 
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OFFICE  OF  SCIENCE  AND  TECH  POLICY 
ATTN:  BHEALY 
ATTN:  COLS  WYMAN 

OFFICE  OF  TECHNOLOGY  ASSESSMENT 
ATTN:  R  WILLIAMSON 

U  S  ARMS  CONTROL  &  DISARMAMENT  AGCY 
ATTN:  B  DOENGES  NWC-DPA 
ATTN:  CDRP  JAMISON 
ATTN:  COL  H  HERTEL 
ATTN:  G  PITMAN 
ATTN:  H  SCHAEFFER 
ATTN:  LTCOL  S  LAWRENCE 
ATTN:  R  GODESKY 
ATTN:  R  HOWES 
ATTN:  R  O’CONNELL  NWC-DPA 

U  S  DEPARTMENT  OF  STATE 
ATTN:  ACORTE 
ATTN:  C  CLEMENT 
ATTN:  COL  M  SEATON 
ATTN:  S  CLEARY 
ATTN:  T  VREBALOVICH 

U  S  GEOLOGICAL  SURVEY 
ATTN:  R DECKER 

U  S  GEOLOGICAL  SURVEY 

ATTN:  E  SHOEMAKER 

U  S  HOUSE  OF  REPRESENTATIVES 
ATTN:  C  BAYER 

ATTN:  COMMITT  ON  SCI  &  TECH  J  DUGAN 

U  S  HOUSE  OF  REPRESENTATIVES 
ATTN:  J  FREIWALD 
ATTN:  MHERBST 

US  DEPARTMENT  AGRICULTURE 
ATTN:  DWARD 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AERO-CHEN  RESEARCH  LABS.  INC 
ATTN:  DB  OLSON 

AERODYNE  RESEARCH,  INC 
ATTN:  C  KOLB 
ATTN:  J  LURIE 

AEROJET  ELECTRO-SYSTEMS  CO 
ATTN:  AFYMAT 
ATTN:  S  HAMILTON 
ATTN:  R  PAN 

AEROSPACE  CORP 

ATTN:  CRICE 
ATTN:  LR  MARTIN 
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AEROSPACE  CORPORATION 
ATTN:  G  LIGHT 

ALLEN  RESEARCH  CORP 
ATTN:  R  ALLEN 

AMERICAN  ASSN  ADVANCEMENT  OF  SCIENCE 
ATTN:  DM  BURNS 

ANALYTIC  SERVICES,  INC  (ANSER) 

ATTN:  RBROFFT 
ATTN:  RELLINGSON 

APPLIED  RESEARCH  CORP 
ATTN:  AENDAL 

ASSN.  DIST.  AMERICAN  SCIENTISTS 
ATTN:  J  HUBBARD 

AT&T  DEFENSIVE  SYSTEMS  STUDIES 
ATTN:  RJANOW 

ATMOSPHERIC  AND  ENVIRONMENTAL  RES 
ATTN:  NS2E 

AUDIO  INTELLIGENCE  DEVICES  INC 
ATTN:  H  BAUM 

AVCO  SYSTEMS  DIVISION 
ATTN:  G  GRANT 

BALL  AEROSPACE  SYSTEMS  DIVISION 
ATTN:  B  CUMMINGS 
ATTN:  C BRADFORD 

BDM  CORP 

ATTN:  D  SHAEFFER 
ATTN:  EL  COFFEY 
ATTN:  J  LEECH 

BERKELEY  RSCH  ASSOCIATES.  INC 
ATTN:  S  BRECHT 

BOEING  AEROSPACE  COMPANY 
ATTN:  NGERONTAKIS 

BOEING  TECHNICAL  &  MANAGEMENT  SVCS,  INC 
ATTN:  GHALL 

C.  L.  CONSULTING  SERVICES 
ATTN:  FFEER 

CALIFORNIA  RESEARCH  &  TECHNOLOGY,  INC 
ATTN:  M  ROSENBLATT 
ATTN:  RGAJ 
ATTN:  S  KRUEGER 

CALSPAN  CORP 

ATTN:  R  MAMBRETTI 
ATTN:  R  MISSERT 


CARNEGIE  CORPORATION  OF  NEW  YORK 
ATTN:  D  ARSENIAN 

CASSIDY  AND  ASSOCIATES 
ATTN:  J JACOBSON 

CHARLES  STARK  DRAPER  LAB,  INC 
ATTN:  ATETEWSKI 

COLORADO  STATE  UNIVERSITY 
ATTN:  D KRUEGER 
ATTN:  W  COTTON 

COMPUTER  SCIENCES  CORP 
ATTN:  G CABLE 

DARTSIDE  CONSULTING 

ATTN:  A  FORESTER 

DELTA  RESEARCH 

ATTN:  L  WEINER 
ATTN:  MRADKE 

DYNAMICS  TECHNOLOGY,  INC 
ATTN:  DHOVE 

ENW  INTERNATIONAL,  LTD 
ATTN:  J  CANE 

EOS  TECHNOLOGIES,  INC 
ATTN:  B  GABBARD 
ATTN:  N JENSEN 
ATTN:  WLELEVIER 

FACTORY  MUTUAL  RESEARCH  CORP 
ATTN:  M  A  DELICHATSIOS 

FEDERATION  OF  AMERICAN  SCIENTISTS 
ATTN:  J  STONE 

GENERAL  ELECTRIC  CO 

ATTN:  RE  SCHMIDT 

GENERAL  ELECTRIC  CO 

ATTN:  H ROBSON 

GENERAL  RESEARCH  CORP 
ATTN:  B  BENNETT 
ATTN:  JBALTES 

HAROLD  ROSENBAUM  ASSOCIATES,  INC 
ATTN:  G  WEBER 

HORIZONS  TECHNOLOGY  INC 
ATTN:  A  EDWARDS 
ATTN:  J  A  MANGO 
ATTN:  J  AMBROSE 

HORIZONS  TECHNOLOGY,  INC 
ATTN:  RWLOWEN 
ATTN:  WTKREISS 
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HUGHES  AIRCRAFT 

ATTN:  EDIVITA 

INFORMATION  SCIENCE.  INC 
ATTN:  W  DUDZIAK 

INSTITUTE  FOR  DEFENSE  ANALYSES 
ATTN:  C  CHANDLER 
ATTN:  E  BAUER 
ATTN:  FALBINI 

JOHNS  HOPKINS  UNIVERSITY 
ATTN:  M  LENEVSKY 
ATTN:  RFRISTROM 
ATTN:  WBERL 

KAMAN  SCIENCES  CORP 
ATTN:  J  RUSH 
ATTN:  J SCRUGGS 

KAMAN  SCIENCES  CORP 
ATTN:  P  GRIFFIN 
ATTN:  P TRACY 

KAMAN  TEMPO 

ATTN:  BGAMBILL 
ATTN:  D  FOXWELL 
ATTN:  DASIAC 
ATTN:  E  MARTIN 
ATTN:  R RUTHERFORD 
ATTN:  R  YOUNG 
ATTN:  SFIFER 
ATTN:  W KNAPP 

KAMAN  TEMPO 

ATTN:  D  ANDERSON 
ATTN:  DASIAC 

LOCKHEED  MISSILES  &  SPACE  CO.  INC 
ATTN:  ATTN  J  HENLEY 
ATTN:  JGLADIS 
ATTN:  J  PEREZ 

LOCKHEED  MISSILES  &  SPACE  CO,  INC 
ATTN:  P  DOLAN 
ATTN:  W  MORAN 

M  I  T  LINCOLN  LAB 

ATTN:  S  WEINER 

MARTIN  MARIETTA  DENVER  AEROSPACE 
ATTN:  D  HAMPTON 

MAXIM  TECHNOLOGIES.  INC 
ATTN:  J MARSHALL 

MCDONNELL  DOUGLAS  CORP 
ATTN:  RC  ANDREWS 
ATTN:  TCRANOR 
ATTN:  TTRANER 


MCDONNELL  DOUGLAS  CORP 
ATTN:  A  MONA 
ATTN:  F  SAGE 
ATTN:  GBATUREVICH 
ATTN:  J  GROSSMAN 
ATTN:  RHALPRIN 
ATTN:  S  JAEGER 
ATTN:  WYUCKER 

MERIDIAN  CORP 

ATTN:  E  DANIELS 
ATTN:  FBAITMAN 

MIDWEST  RESEARCH  INSTITUTE 
ATTN:  JS  KINSEY 

MISSION  RESEARCH  CORP 

ATTN:  R  ARMSTRONG 

MISSION  RESEARCH  CORP 
ATTN:  C  LONGMIRE 
ATTN:  D ARCHER 
ATTN:  DKNEPP 
ATTN:  DSOWLE 
ATTN:  FFAJEN 
ATTN:  KRCOSNER 
ATTN:  M  SCHEIBE 
ATTN:  RBIGONI 
ATTN:  R  CHRISTIAN 
ATTN:  R  GOLDFLAM 
ATTN:  R  HENDRICK 
ATTN:  TOLD 
ATTN:  W  WHITE 

MITRE  CORPORATION 

ATTN:  J  SAWYER 

MRJ  INC 

ATTN:  D  FREIWALD 

NATIONAL  ADVISORY  COMMITTEE 
ATTN:  JALMAZAN 
ATTN:  J  BISHOP 

NATIONAL  INST.  FOR  PUBLIC  POLICY 
ATTN:  K  PAYNE 

NICHOLS  RESEARCH  CORP,  INC 
ATTN:  H  SMITH 
ATTN:  J  SMITH 
ATTN:  M  FRASER 
ATTN:  RBYRN 

NORTHROP  SERVICES  INC 
ATTN:  T OVERTON 

ORLANDO  TECHNOLOGY  INC 
ATTN:  R  SZCZEPANSKI 
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PACIFIC-SIERRA  RESEARCH  CORP 
ATTN:  GANNO 

ATTN:  H  BRODE.  CHAIRMAN  SAGE 
ATTN:  MOORE 
ATTN:  R  SMALL 

PALOMAR  CORP 

ATTN:  B GARRETT 
ATTN:  CFELDBAUM 

PHOTOMETRICS.  INC 

ATTN:  ILKOFSKY 

PHOTON  RESEARCH  ASSOCIATES 
ATTN:  JMYER 

PHYSICAL  RESEARCH  CORP 
ATTN:  ACECERE 

PHYSICAL  RESEARCH  INC 
ATTN:  HFITZ 

PHYSICAL  RESEARCH  INC 
ATTN:  D  MATUSKA 

PHYSICAL  RESEARCH  INC 

ATTN:  AWARSHAWSKY 
ATTN:  J  WANG 
ATTN:  WSHIH 

PHYSICAL  RESEARCH  INC 
ATTN:  RJORDANO 

PHYSICAL  RESEARCH.  INC 
ATTN:  DWESTPHAL 
ATTN:  DWHITENER 
ATTN:  H  WHEELER 
ATTN:  R  BUFF 
ATTN:  R  DELIBERIS 
ATTN:  T  STEPHENS 
ATTN:  W  C  BLACKWELL 

PHYSICAL  RESEARCH.  INC 
ATTN:  G  HARNEY 
ATTN:  J  DEVORE 
ATTN:  J  THOMPSON 
ATTN:  R  STOECKLY 
ATTN:  W  SCHLEUTER 

PHYSICAL  RESEARCH.  INC 
ATTN:  HSUGIUCHI 

POLYTECHNIC  OF  NEW  YORK 
ATTN:  BJ  BULKIN 
ATTN:  GTESORO 

PRINCETON  UNIVERSITY 

ATTN:  JMAHLMAN 

QUADRI  CORP 

ATTN:  H  BURNSWORTH 


R  &  0  ASSOCIATES 

ATTN:  AKUHL 
ATTN:  F  GILMORE 
ATTN:  G  JONES 
ATTN:  J SANBORN 
ATTN:  RTURCO 

R  &  D  ASSOCIATES 

ATTN:  8  YOON 

R  J  EDWARDS  INC 

ATTN:  R  SEITZ 

RADIATION  RESEARCH  ASSOCIATES.  INC 
ATTN:  B  CAMPBELL 
ATTN.  M  WELLS 

RAND  CORP 

ATTN:  GL  DONOHUE 
ATTN:  P  ROMERO 

RAND  CORP 

ATTN:  JGERTLER 

ROCKWELL  INTERNATIONAL  CORP 
ATTN:  S  I  MARCUS 

ROCKWELL  INTERNATIONAL  CORP 
ATTN:  J  KELLEY 

S-CUBED 

ATTN:  B  FREEMAN 
ATTN:  KDPYATT.JR 
ATTN:  R  LAFRENZ 

S-CUBED 

ATTN:  C NEEDHAM 
ATTN:  SHIKIDA 
ATTN:  T  CARNEY 

SCIENCE  APPLICATIONS  INC 
ATTN:  REDELMAN 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  CHILL 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  D  HAMLIN 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  B  MORTON 
ATTN:  B  SCOTT 
ATTN:  D  SACHS 
ATTN:  GT  PHILLIPS 
ATTN.  J  BENGSTOM 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  D  BACON 
ATTN:  DRLGOURE 
ATTN:  FGIESSLER 
ATTN:  J  COCKAYNE 


DASIACTN- 86-29-V3  (DL  CONTINUED) 


ATTN:  J  SHANNON 
ATTN:  J STUART 
ATTN:  MSHARFF 
ATTN:  WLAYSON 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  J  SONTOWSKI 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  T  HARRIS 

SCIENTIFIC  RESEARCH  ASSOC,  INC 
ATTN:  B  WEINBERG 

SPARTA  INC 

ATTN:  R HARPER 

SRI  INTERNATIONAL 

ATTN:  CWITHAM 
ATTN:  D  GOLDEN 
ATTN:  D  MACDONALD 
ATTN:  D  ROBERTS 
ATTN:  EUTHE 
ATTN:  GABRAHAMSON 
ATTN:  J  BACKOVSKY 
ATTN:  WCHESNUT 
ATTN:  W JOHNSON 

SRI  INTERNATIONAL 

ATTN:  R  BRAMHALL 
ATTN:  R  WOOLFOLK 
ATTN:  W  VAIL 

STAN  MARTIN  ASSOCIATES 
ATTN:  SB  MARTIN 

STANTON  CONSULTING 

ATTN:  M STANTON 

SWETL,  INC 

ATTN:  TY  PALMER 

SYSTEM  PLANNING  CORP 
ATTN:  J  SCOUR  AS 
ATTN:  MBIENVENU 
ATTN:  R  SCHEERBAUM 

SYSTEMS  AND  APPLIED  SCIENCES  CORP 
ATTN:  M  KAPLAN 

TECHNOLOGY  INTERNATIONAL  CORP 
ATTN:  WBOQUIST 

TELEDYNE  BROWN  ENGINEERING 
ATTN:  D  ORMOND 
ATTN:  F LEOPARD 
ATTN:  J  FORD 

TELEDYNE  BROWN  ENGINEERING 
ATTN:  DGUICE 


TEXAS  ENGR  EXPERIMENT  STATION 
ATTN:  WH  MARLOW 

TOYON  RESEARCH  CORP 
ATTN:  CTRUAX 
ATTN:  J  GARBARINO 
ATTN:  J  ISE 

TRW 

ATTN:  H  BURNSWORTH 
ATTN:  JBELING 

TRW  ELECTRONICS  &  DEFENSE  SECTOR 
ATTN:  FFENDELL 
ATTN:  G  KIRCHNER 
ATTN:  G  MROZ 
ATTN:  H  CROWDER 
ATTN:  J  FEDELE 
ATTN:  M  BRONSTEIN 
ATTN:  R  BACHARACH 
ATTN:  S  FINK 
ATTN:  T  NGUYEN 

TRW  ELECTRONICS  &  DEFENSE  SECTOR 
ATTN:  M  HAAS 

VISIDYNE,  INC 

ATTN:  H  SMITH 
ATTN:  J CARPENTER 

WASHINGTON,  UNIVERSITY  OF 
ATTN:  J  I  KATZ 

FOREIGN 

AERE  ENVIRONMENTAL  AND  MEDICAL  SC 
ATTN:  SPENKETT 

ATOMIC  WEAPONS  RESEARCH  ESTABLISHMENT 
ATTN:  PF  A  RICHARDS 

ATOMIC  WEAPONS  RESEARCH  ESTABLISHMENT 
ATTN:  DL  JONES 
ATTN:  DM  MOODY 

AUSTRALIA  EMBASSY 

ATTN:  DR  LOUGH 
ATTN:  MAJ  GEN  H  J  COATES 
ATTN:  P PROSSER 

BRITISH  DEFENCE  STAFF 
ATTN:  C  FENWICK 
ATTN:  JCRANIDGE 
ATTN:  J  EDMONDS 
ATTN:  M  NORTON 
ATTN:  P  WEST 

CANADIAN  FORESTRY  SERVICE 
ATTN:  B  STOCKS 
ATTN:  TLYNHAM 
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CSIRO 

ATTN:  I  GALBALLY 

CSIRO:  ATMOSPHERIC  RESEARCH 
ATTN:  APITTOCK 

EMBASSY  OF  BELGIUM 

ATTN:  LARNOULD 

ISRAEL  EMBASSY 

ATTN:  N  BELKINO 

MAX-PLANCK  INSTITUTE  FOR  CHEMISTRY 
ATTN:  P  J  CRUTZEN 

MINISTRY  OF  DEFENCE 
ATTN:  R  RIDLEY 

NATIONAL  DEFENCE  HEADQUARTERS 
ATTN:  HAROBITALLE 

TRINITY  COLLEGE 

ATTN:  F  HARE 

DIRECTORY  OF  OTHER 

ATMOS.  SCIENCES 

ATTN:  GSISCOE 

BROWN  UNIVERSITY 

ATTN:  RK  MATTHEWS 

BUCKNELL  UNIVERSITY 

ATTN:  0  ANDERSON 

CALIFORNIA.  UNIVERSITY 

ATTN:  R  WILLIAMSON 

CALIFORNIA.  UNIVERSITY  OF 

ATTN:  L  BADASH/DEPT  OF  HISTORY 

COLORADO.  UNIVERSITY  LIBRARIES 
ATTN:  JBIRKS 
ATTN:  R  SCHNELL 

DREXEL  UNUVERSITY 

ATTN:  J  FRIEND 

DUKE  UNIVERSITY 

ATTN:  FDELUCIA 

GEORGE  MASON  UNIVERSITY 
ATTN:  PROFS  SINGER 
ATTN:  R  EHRLICH 

GEORGE  WASHINGTON  UNIVERSITY 
ATTN:  R  GOULARD 

GEORGIA  INST  OF  TECH 

ATTN:  E  PATTERSON 


HARVARD  COLLEGE  LIBRARY 
ATTN:  W  PRESS 

HARVARD  UNIVERSITY 

ATTN:  G  CARRIER 

HARVARD  UNIVERSITY 

ATTN:  DEARDLEY 

IOWA,  UNIVERSITY  OF 

ATTN:  HISTORY  DEPT/3  PYNE 

MARYLAND  UNIVERSITY  OF 

ATTN:  A  ROBOCK  DEPT  METEOROLOGY 
ATTN:  A  VOGELMANN  DEPT  METEOROLOGY 
ATTN:  R  ELLINGSON  DEPT  METEOROLOGY 

MIAMI  LIBRARY  UNIVERSITY  OF 
ATTN:  C  CONVEY 

MIAMI  UNIV  LIBRARY 

ATTN:  J  PROSPERO  ATMOS  SC 

NEW  YORK  STATE  UNIVERSITY  OF 
ATTN:  R  CESS 

OAK  RIDGE  ASSOCIATED  UNIVERSITIES 
ATTN:  C  WHITTLE 

PENNSYLVANIA  STATE  UNIVERSITY 
ATTN:  DWESTPHAL 

SOUTH  DAKOTA  SCH  OF  MINES  &  TECH  LIB 
ATTN:  H  ORVILLE 

TENNESSEE.  UNIVERSITY  OF 
ATTN:  K  FOX 

UNIVERSITY  OF  SOUTH  FLORIDA 
ATTN:  SYING 

UNIVERSITY  OF  WASHINGTON 
ATTN:  CLEOVY 
ATTN:  LRAOKE 
ATTN:  P  HOBBS 

VIRGINIA  POLYTECHNIC  INST  LIB 
ATTN:  MNADLER 

WASHINGTON  STATE  UNIVERSITY 
ATTN:  DR  A  CLARK 

WISCONSIN  UNIVERSITY  OF 
ATTN:  P  WANG 
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